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Abstract 

We study the addition of localised impurities to U(N) Super symmetric Yang-Mills theories 
in (p+ l)-dimensions by using the gauge/gravity correspondence. From the gravity side, the 
impurities are introduced by considering probe D(8 — p)-branes extending along the time and 
radial directions and wrapping an (7 — p)-dimensional submanifold of the internal (8 — p)- 
sphere, so that the degrees of freedom are point-like from the gauge theory perspective. We 
analyse both the configuration in which the branes generate straight flux tubes - correspond- 
ing to actual single impurities - and the one in which connected flux tubes are created - 
corresponding to dimers. We discuss the thermodynamics of both the configurations and the 
related phase transition. In particular, the specific heat of the straight flux-tube configura- 
tion is negative for p < 3, while it is never the case for the connected one. We study the 
stability of the system by looking at the impurity fluctuations. Finally, we characterise the 
theory by computing one- and two-point correlators of the gauge theory operators dual to 
the impurity fluctuations. Because of the underlying generalised conformal structure, such 
correlators can be expressed in terms of an effective coupling constant (which runs because 
of its dimensionality) and a generalised conformal dimension. 
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1 Introduction 



The gauge/gravity correspondence [1-4] provides a powerful handle on the strongly coupled 
regime of gauge theories. As originally formulated [1], it conjectures a relation between the 
AdSs x S 5 "near-horizon" geometry generated by a stack of N coincident D3-branes and the 
four-dimensional N = 4 SU(N) Supersymmetric Yang-Mills theory living on the boundary 
of AdS§. Such an equivalence between supergravity on AdS$ x S 5 on one side and Af = 4 
SYM theory on the other, can be made precise by identifying the supergravity partition 
function with the generating function for the gauge theory correlators, with the boundary 
value of the bulk modes acting as a source of the corresponding gauge theory operator [3]. 

The correspondence can be extended to the case of generic Dp-branes (p ^ 3), whose 
"near-horizon" geometry describes the adjoint degrees of freedom of U(N) Supersymmetric 
Yang-Mills theory in (p+ l)-dimensions [5]. One main difference with the p = 3 case is that, 
while the latter is characterised by a dimensionless coupling constant, the coupling constant 
turns out to be dimensionful for p / 3 and the effective coupling constant runs with the 
energy scale, making the theory non-conformal. However, there exists a frame [6] - named 
dual frame - in which the background induced by the stack of N Dp-branes is conformally 
AdSp+i x S s ~ p with the conformal factor depending on a non-trivial dilaton [7-9]. This frame 
is further characterised by a manifest generalised conformal symmetry [10-12]. Specifically, 
the generalised conformal structure emerges by allowing the string/Yang-Mills coupling to 
transform as a background field under conformal transformation, providing, on the other 
hand, diffeomorphism and trace Ward identities. Moreover, in this frame the radial direction 
plays the role of the energy scale of the dual gauge-theory [9,13]. The holographic RG flow is, 
in any case, trivial since the theory flows just because of the dimensionality of the coupling 
constant. The dual frame has been crucial to make holography for Dp-branes precise and 
to extend the holographic renormalisation procedure 1 [15-17]. In particular, it has been 
observed [17] that the (p + 2)-effective supergravity action, obtained by the Kaluza-Klein 
reduction on the (8 — p)-dimensional sphere, can be recovered by dimensional reduction of 
the theory on pure asymptotically AdS2 a +i on a torus, where a is a parameter related to the 
power of the radial direction in the dilaton and, generically, takes fractional values for generic 
values of p. Such a parameter a can be considered as an integer and, after the reduction on the 
torus, can be analytically continued to take its actual p-dependent value. One can therefore 
map the problem to a pure AdSi<j+\ theory and then obtain the (p + 2)-dimensional answer 
through Kaluza-Klein reduction. This way to rephrase the problem allows to drastically 
simplify the computation of the counterterms needed for holographic renormalisation and of 
quantities such as, for example, the correlators of the stress-energy tensor. 
1 For more details about holographic renormalisation, see [14]. 
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As we mentioned earlier, the backgrounds generated by a stack of Dp-branes describe 
the adjoint degrees of freedom of the dual gauge theory. The correspondence can however 
be generalised by inserting extra degrees of freedom in the theory. In particular, one can 
add a certain number of branes as probes, so that they do not backreact on the geometry, 
introducing a fundamental hyper-multiplet and partially or completely breaking the initial 
super symmetries [18]. 

In this paper we are concerned with the introduction of D(8 — p)-branes in the conformally 
AdSp+i x S 8 ~ p background generated by a stack of Dp-branes, in such a way that they wrap 
an 5 7 ~ p C S 8 ~ p and the induced metric on their world-volume turns out to be conformally 
AdS2 x S 7 ~ p . Such configurations were analysed in [19, 20] and they describe bound states 
of fundamental strings stretching along the radial direction. From the point of view of the 
dual SYM theory, this is equivalent to introducing point like objects acting as impurities in 
an interacting (trivially) non-conformal field theory. The conformal case p = 3 has been 
recently studied to holographically realise Kondo models [21-25], and such idea has been 
extended to the case of ABJM theory for which the impurities have been introduced via D6- 
branes extended along AdS2 C AdS± and wrapping a squashed T 1, space inside the internal 
CP 3 [26]. 

We mainly study two classes of configurations: one in which a straight flux tube is formed 
and the other in which two flux tubes ending on two different points at the boundary are 
connected in the bulk. In order to establish the stability of the brane configurations, we 
analyse the fluctuations of the probe branes at zero temperature. After decoupling the 
modes, we can characterise the operators dual to such fluctuations through the effective 
coupling constant and the generalised scaling dimension. In particular, we can read them 
off by computing one-point and two-point correlation functions in the coordinate space. The 
calculation is made straightforward by observing that such fluctuations satisfy the equation 
of motion of a free massive scalar propagating in a higher dimensional AdSq+i space, where 
the dimension q needs to be analytically continued to a (generically) fractional value 2 . The 
dual operators are irrelevant, so we can make sense of the correlation functions by using the 
observation made in [29] stating that it is possible to holographically renormalise the bulk 
fields in a perturbative way up to some fixed order n. Given that we are interested into one- 
point and two-point correlators, it is enough to holographically renormalise the fluctuation 
action up to quadratic order, where it can always be written as the action of a free massive 
scalar in AdS q +\. The results we find for the generalised scaling dimensions coincide with 
the known conformal dimensions for such operators in the case p = 3. 

2 This has already been observed in the case of probe flavour branes [27] in these non-conformal backgrounds, 
for which mainly the embedding functions have been considered. For a discussion about the holographic 
renormalisation of probe flavour branes in the conformal case, see [28]. 
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We further analyse the basic thermodynamic properties of the impurities, such as their 
free energy and entropy. Interestingly, we find that, in the class of systems of interest, the 
impurity entropy is generically non-analytic in the filling fraction, except for the D4/D4- 
system for which it is possible to obtain a simple closed expression. This is actually the only 
system in which the impurity specific heat turns out to be positive, while it is zero for p = 3 
and it is negative for p < 3. 

The paper is organised as follows. In Section 2 we review the basic features of the 
Kondo models. In particular, we discuss the effect of a single impurity-spin in a gas of free 
fermions in (3 + l)-dimensions with one or k identical channels and with an SU(2) or SU(N) 
spin-symmetry group. We also review the case in which the impurity spin is inserted in 
(1 + l)-dimensional systems, in which the interactions cannot be generally neglected. In Sec- 
tion 3 we discuss the holographic setup for the ambient non-conformal theory, emphasising 
the existence of a frame in which the bulk is conformally AdS, and we also introduce the 
impurities as probe D(8 — p)-branes wrapping a conformally AdS2 x S 7 ~ p submanifold of the 
bulk geometry. We also comment on the underlying generalised conformal structure in the 
systems we are considering. In Sections 4 and 5 we discuss the two classes of configurations 
of interests, respectively the one in which the probe branes generate straight flux tubes and 
the one in which two flux tubes ending on two different points at the boundary are connected 
in the bulk. In particular, Section 4 is dedicated to the study of the thermodynamics of the 
straight-flux tube/impurity by computing explicitly the free energy, the entropy, the internal 
energy, the specific heat and the susceptibility as functions of the temperature and of the 
filling fraction. Section 5 is instead dedicated to a preliminary study of the hanging flux tube 
configurations, for which the free energy and the one-point function of the operator dual to 
the embedding function are computed. In Section 6 we study the thermodynamics of such 
connected configurations. We again compute the free energy, entropy, internal energy, and 
specific heat and we study the competition between the dimer configurations and the discon- 
nected impurities. We identify a first order phase transition and the related temperature at 
which it occurs. Section 7 is devoted to the detailed analysis of the fluctuations for the probe 
branes. We identify two decoupled channels and we compute the one-point and two-point 
correlators of the operators dual to such fluctuations, emphasising their generalised conformal 
structure. Finally, Section 8 contains our conclusion. 

2 Generalities of the Kondo Model 

The Kondo model [30-33] is one of the examples of quantum impurity problems [33] , where 
by impurity one refers to point-like degrees of freedom inserted in and interacting with a 
(generally non-interacting) gas. This class of systems is characterised by two main features: 
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the existence of gapless excitations far away from the impurity and the localisation at one 
point in position space of the interaction with the impurity. Furthermore, even if the impuri- 
ties are introduced in (generally) (3 + l)-dimensional theories, such systems can be mapped 
to (1 + l)-dimensional one - except for the case in which the bulk theory is already (1 + 1)- 
dimensional -, and the impurities appear located at the origin of the one-dimensional position 
space. 

Typically, one talks about Kondo problem when a free fermion gas is coupled to an 
impurity with a spin degree of freedom. The interesting feature of such a system is that 
the resistivity shows a minimum as a function of the temperature: when the temperature 
approaches the so-called Kondo temperature, the interaction of the bulk electrons with the 
impurity starts to compensate the effect of the interaction with the lattice phonons (which 
leads to a decrease of the resistivity as the temperature is lowered) so that, as the temperature 
becomes lower than the Kondo temperature, the effect of the impurity dominates. 

2.1 Single-Channel Kondo Model 

The Hamiltonian density describing the Kondo model with just one impurity spin is 



v/ 2 i , ,a 

W (s +i) = K 



2m 



llfa + JSM&Wl—fa-t, (2.1) 



where ip a is the bulk electron field with spin- index a, e F is the Fermi energy, a are the Pauli 
matrices, S is the spin-1/2 impurity operator and J is the coupling for a Heisenberg-type 
exchange interaction between the impurity spin and the electron spin density. Moreover, a 
sum over the spin-index a is understood. 

The coupling constant J is dimensionful. Its dimensionless version can be defined as 
A = f Jp, p being the density of states per unit of energy, per unity volume and per spin. In 
the case of a free fermion gas p = mk F /ir 2 , with k F being the Fermi momentum. 

The delta function in the coupling selects just the s-wave to interact with impurity. This 
allows the (3+l)-dimensional model described by the Hamiltonian density (2.1) to be reduced 
to an effective (1 + l)-dimensional model defined on half-line with the impurity located at 
the origin: 

n (1 + 1) =^ a i^ a + X6(x)^ a —^ a -~i. (2.2) 

This is the Hamiltonian density for a massless Dirac fermion in (1 + l)-dimensions and 
therefore the theory is conformal, with a boundary. A characteristic of this system is that 
the renormalised coupling constant X(E) increases as the energy scale E is lowered, so that 
at low enough energies perturbation theory breaks down. The /3-function turns out to be 

-A 2 + -.., =► A(A) ~ - - , A ° A - , (2.3) 



din A " ' w l-A ln(A /A) ' 
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A being a scale, and A is the scale at which the Kondo coupling constant acquires the bare 
value A . In the case in which the bare coupling A is negative (ferromagnetic behaviour), the 
coupling A decreases as the energy scale decreases, being therefore well-behaved. If instead 
A is positive, there is a value for the energy scale at which perturbation theory breaks down 
and it is represented by the so-called Kondo temperature T K 

T K ~ Aoe- 1 /*". (2.4) 

In this case, the Kondo coupling constant can be thought to renormalise to infinity [34-36]. 
This can be understood by looking at the lattice version of the (l+l)-dimensional Hamiltonian 
density (2.2) 

U\? + 1} = -t + h -c. ] + J 4 y^o • 3, (2.5) 

i 

where the strong coupling regime is at J S> t. For t = 0, the electron configuration in 
a general site is arbitrary, except at the origin where the electron form a singlet with the 
impurity. For relatively small t compared with the Kondo coupling A, the electrons in a 
general site are in a Bloch state whose single-particle wave-function vanishes at the origin in 
order to preserve the singlet-condition. 

The interesting feature of this system is that, in the even-parity sector, the wave- functions 
at zero and infinity Kondo couplings differ from each other, while in the odd-parity sector 
they do not. The strong coupling fixed point is the same as the weak coupling one, with the 
difference that the impurity is screened and substituted by the boundary condition ip(0) = 0. 
For finite/small Kondo coupling this is still true but only at low energies and long distances. 
The boundary condition is a fixed point and this is a feature of all quantum impurity systems. 

The flow to the low energy fixed point is controlled by the leading irrelevant operator, 
which is constructed out of the fermion fields and it is 5C/(2)-invariant. Such a symmetry al- 
lows for two dimension-2 operators, (ip^ a (0)ip a (0)) 2 and (^(ip^ a (0)~c? a ^ip a (0))/2 S J . However, 
the first one is suppressed given that its coefficient turns out to be of order 1 /A , while the 
coefficient of the latter operator of order 1/T K , with A S> T K . 

Using perturbation theory in 1/T K , one finds that, at low temperature, the impurity 
susceptibility - defined as the difference between susceptibility with and without impurity - 
turns out to behave as a constant. Furthermore, at high temperature and in the scaling limit 
of small bare Kondo coupling, it is proportional to the inverse of the temperature. Actually, 
the impurity susceptibility can be written as 



imp 1 f T 



x i»P( T ) _^ (4T K )-\ as T -> , 

(2.6) 

X imp (T) — > (AT)- 1 , for T > T K . 
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2.2 Mult i- Channel Kondo Model 



The single-impurity model discussed in the previous section can be generalised by considering 
k conduction electron channels interacting with the impurity spin operator S [37, 38] 



+ 



(2.7) 



where a is a spin-index (for the moment considered to run on two possible spin-states), the 
sum over the index j is understood and j = l,...,k. In this Hamiltonian density, the 
channels are assumed to be identical, preserving a SU (fc)-symmetry. 
The /3-function turns out to be 



dlnA " " A + 2 + " 



(2.8) 



As in the single-channel case, the Kondo coupling A renormalises to infinity for an anti- 
ferromagnetic bare coupling A , as long as S > k/2. More precisely, one electron per channel 
is expected to go into a symmetric state near the origin forming a total spin k/2 and the 
ground state has size \S — k/2\. The low temperature features in the under-screened case 
S > k/2 are similar to the ones of the single-channel Kondo model, where the strong coupling 
fixed point is stable. The same occurs in the exactly screened case S = k/2. In the over- 
screened case S < k/2, the RG-flow leads to a non-trivial finite coupling fixed point, for 
which the free energy is non-analytic as function of temperature and magnetic field. Let us 
look more in detail to the thermodynamic properties. In particular, the impurity specific 
heat and the impurity magnetic susceptibility in the over-screened case turn out to be 



C imp {T) oc 



rp4/{k + 2) 



^ln^ 

T K 111 T K ■ 



k ^ 2 



k 



X imp (T) oc 



fe— 2 



In 2k 

111 rp , 



k ^ 2 



(2.J 



In the under-screened case, the impurity magnetic susceptibility diverges as the temperature 
T is lowered, leading to a non-Fermi-liquid behaviour. In the large-A: case, instead, one has 



C imp {T) oc 1 + -lnT + O (AT 2 ) 



T 



1 + - InT + (k~ 2 ) 

Ki 



(2.10) 



The magnetic susceptibility diverges as T~ l as the temperature is lowered, while a low- 
temperature logarithmic divergence in the impurity specific heat appears as a l//c-effect. 

Let us now consider the case in which the spin-symmetry group is SU(N), so that the 
spin-index a in the Hamiltonian density (2.7) run from 1 to N [39]. The bulk fermions trans- 
form under the fundamental representation of SU(N), while the impurity-spin transforms 
under an anti-symmetric representation with Q indices. As in the SU (2) multichannel model 
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discussed above, also its SU(N) generalisation presents a fixed point at intermediate Kondo 
coupling. In this fixed point, the local impurity spin two-point function has the following 
zero-temperature behaviour 

< S (^(0)> ~ A imp = j^L., (2.11) 

and the impurity susceptibility behaves as the local susceptibility (which is obtained inte- 
grating the above correlation function) 

k-N 

T ~, k > N, 



X imp (T) ~ X loc (T) ~ { 



lnT" 1 , k = N, (2.12) 



const., k < N 



Notice that the expression (2.12) is the natural generalisation to the SU(N)-spin symmetry 
of Eq. (2.9), where the spin symmetry is SU(2) - the impurity susceptibility acquires a 
logarithmic behaviour when the number of channels is equal to the rank of the spin symmetry 
group. 

It is interesting to consider the l&rge-N limit of this class of systems. There are two 
ways to take such a limit: one can take just the rank of the spin-symmetry group to be 
large (with k <C N) or taking it to be large and keeping the ratio 7 = f k/N between the 
number of channels and the rank of the spin-symmetry group to be fixed. In the latter case, 
the expression for the scaling dimension of the local impurity spin and the behaviour of the 
impurity susceptibility with the temperature do not change - one can just conveniently re- 
write them in terms of the parameter 7 as A imp = (1 + 7) _1 and x imp ~ T~( 7_1 )/( 7+1 ) for 
7 > 1. In case k is instead not taken to be large (and of the same order of N), the large- N 
limit of the local impurity spin scaling dimension and of the impurity susceptibility become 

A imp = 1 - ^ + O ((k/N) 2 ) , X imp (T) ~ const., k <C N. (2.13) 



N 

Finally, it is interesting to write down the explicit expression for the impurity entropy 

Q 



S imp = In M — , (2.14 

sin [tt v/ ( k + TV )J 



which vanishes in the single-channel case. 



2.3 Impurities in Luttinger Liquids 

The impurity systems reviewed in the previous two sections are characterised by the fact that 
the bulk fermions are non-interacting and the localisation of the impurity at a point allows to 
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reduce the original (3 + l)-dimensional problem to a (1 + l)-dimensional one. However, one 
can think to introduce an impurity spin directly in a (1 + l)-dimensional system. In this case, 
the bulk degrees of freedom cannot always be considered as non-interacting and, moreover, 
they cannot behave neither as a Bose liquid nor as a Fermi liquid. Rather, they are thought 
to be described by Tomonaga-Luttinger liquids [40-45]. 

This class of systems is characterised by the absence of quasi-particle excitations and the 
presence of plasmons and spin density waves, which are independent of each other. Further- 
more, while for Bose and Fermi liquids the specific heat at low temperature scales as T d , 
d > 1 being the number of spacial dimensions, and T respectively, in (1 + l)-dimensions it 
can scale either as ~ T Q ( A(L) ) or linearly with the temperature, depending on whether the 
system is interacting or not - a(A (L) ) is a function of the Luttinger parameter A (L) and its 
form depends on the attractive or repulsive nature of the interaction. 

The Hamiltonian density for a single impurity in a Luttinger liquid can be written as 

n Lut = -v F ^ [(v>jV>j+i + h.c. ) + utij] + ?C k + ^ imp , (2-15) 

3 

where the first term is the Hubbard model describing the bulk degrees of freedom, with fij 
being the total- number operator at site j, v F = 2at sin (ak F ), U and t are the usual Hubbard 
parameters, %!° n " lk is the bulk interaction term, and T~L luip is the bulk-impurity interaction. 
Several bulk interactions are possible. Examples are the operators 0\ = A (L) J L J R , O2 = 

R , as well as spin anisotropic interactions 
of zero conformal spin. In the case 7i^ k ~ Oi, through a field redefinition, the Luttinger- 
liquid Hamiltonian density acquires the non-interacting form, at the price that the scaling 
dimensions of various operators change. If 'H| 3 n " lk ~ O2 is considered, its effect is not-negligible 
just in the half-filling case and produces a charge gap for values of the Luttinger parameter 
for which the interaction is repulsive. 

Let us now turn to the bulk-impurity term of T~L Lut , which can be taken as 

^ imp = A fc ^L^f, (2.16) 

where the indices k, I run over L,R, which indicates left- and right-movers, while the indices 
a, (3 run over the spin values - in the case of SU{2) spin-symmetry, they can take just two 
values. For A LL = A RR and A LR = A RL one gets respectively the amplitudes for forward and 
backward scattering of the bulk degrees of freedom with the impurity. The Kondo interaction 
is obtained when all these couplings are equal. 

In a model with such a bulk-impurity interaction with the bulk-interaction switched off, 
the Hamiltonian density can be mapped to the one of a two-channel model with the impurity 
coupled to the bulk fermions just in one of the channels. Like in the (3+l)-dimensional 
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Kondo problem, the model renormalises to a local Fermi liquid and the response function 
scales analytically with the temperature [37]. 

When a bulk interaction such as O2 is switched on, any impurity interaction is expected to 
be substituted by a renormalised boundary condition on the critical bulk theory. The theory 
turns out to be characterised by the impurity specific heat and the impurity susceptibility at 
low temperature given by [46] 



C imp (T) = ci 



1 

1 



~i 2 1 1 

' (L) x 



T a \ xK ') + C2 T, x imp (T) = c 3 T°, (2.17) 



witha(A (L) ) = (l + 2A (L Vt>p)" 1/2 . 
3 The Holographic Set-up 

Let us start with a brief review of the holographic realisation of the "ambient" theory where 
the impurities will be later introduced. Consider the background generated by a stack of N 
Dp-branes in the string frame: 



/ r-pX-Va / 7—p\ 1 / 2 

ds 2 w = 1 + ^f- rt^dx" + 1 + \-\ ds 2 T , (3.1) 



where /x, v = 0, . . . ,p, r}av is the flat Minkowski metric in p + 1 dimensions, r is the radial 
coordinate of the transverse space (with ds^ being its line element) and the constant r p is 
defined through 



r 7-p def 
V 



(2V^r P T (IZ*) 9s N(a'f- p)/2 = d p g s N(a'f'^ 2 , (3.2) 



where, in the last step, we have introduced the numerical constant d p , defined as: 



d p = (2V^f- p r( 7 -^-) 



The decoupling limit 



g s -> 0, a' ->■ 0, U = —, = fixed, gl M N = fixed, (3.4) 

a' 

where the coupling constant g Y M is dimensionful and defined by 

gl M d = c 2(2vrr 2 5s (a') (P " 3)/2 5 (3-5) 



corresponds to the "near-horizon" geometry for Dp-branes. In this limit the metric becomes 

lx N = 

\ (7-p)/2 frj x (7-p)/2 ) /o 



ds\ Q = gMNdx M dx N 



a 



U p 



\ \<-V)l* /TT \ K'-P)/^ 
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In (3.6) d£l\_ p is the line element of a S 8 p sphere and we have introduced the constant U p , 
which is defined as: 

r/7-P 4i f d P 2 N (r , 7) 

p ~ 2(2tt)p-2 5ym ' 1 ' 

The Dp-brane background is also endowed with a non-trivial dilaton <f> and a RR (p+ l)-form 

potential C^ +1 ^, which are given by: 

, / 17 \ (7 " p)(p - 3)/4 „(,+!) _ 2(2vrr 2 (aQ^ 1 )/ 2 iV / 17 \ 7 ~ p 

As it has just been mentioned, for p / 3 the coupling constant is dimensionful and, therefore, 
the effective coupling runs with energy scale as: 

g 2 B = g 2 YM NU p ~ 3 . (3.9) 

It turns out that the background metric (3.6) is conformal to an AdS p+ 2 X S 8 ~ p space for 
p / 5 [6,9, 13]. This feature can be seen explicitly by redefining the radial coordinate as 
follows 

and rewriting the line element (3.6) as 

7 2/(7-p) 



2 

dsl = (Ne*y {7 ~ P) B p ds 2 10 , B p d = f a'^j— , (3.11) 

Up 



so that the line element ds 2 describes an AdSp+2 x S 8 p geometry in the Poincare patch 

du 2 

ds 10 = g M Ndx M dx N = u ri llu dx IJ 'dx' / H + u 2 dQ%_ v . (3-12) 

u 

In such a frame, the radial coordinate u plays the role of the energy scale of the boundary 
theory and its rescalings are just the dilatations in the boundary theory [13] . One can also use 
the Fefferman-Graham coordinates by introducing a new radial coordinate p as the inverse 
of the current one, namely 

p = u~ l . (3.13) 
In terms of p the metric ds 2 Q of (3.12) can be written as 

,~2 rj^dx^dx" + dp 2 2 2 

ds 10 = -J hu p dn 8 _ p , (3.14) 

with dilaton and (8 — p)-form field strength F^ 8 ~ p ^ = *dC <yP+l ^ which become 

C " = 2$f*N {^ P P )" L&fl ' Fi8 ~ P) = V~v) dpN ^ ™(S 8 ~ P ) = dCi7 - p) > 

(3.15) 
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Figure 1: The points of the S^ 8 p ^ sphere with the same latitude 9, measured from one of its 
poles, define a S^ 7 ~ p ^ sphere. 



where Vol (S q ) denotes the volume form of a S q round sphere. In order to write explicitly 
the form of the (7 — p)-form potential C^ 7 ~ p \ let us represent the line element of the S^ s ~ p ^ 
sphere in terms of a polar angle 9 and the line element of a S^ 7 ~ p ^ sphere, namely 

dnl_ p = de 2 + sin 2 9 dn^ p , (3.16) 

with 9 taking values in the range < 9 < ir (see figure 1). Then, the (7 — p)-form potential 
can be taken to be 

C 7 _ p = -Nd p {a'f^' 2 C p {9) Vol(S 7 ~ p ) , (3.17) 
where C p (9) is the function uniquely defined by the conditions: 

d £M = _ (7 -p) s in 7 -^, C p (0)=0. (3.18) 

In the case of black hole geometries, the metric (3.14) becomes 

dsl = ~ hpip)d !" + — + K( P )r ld 4 + <dnj_ p , (3.19) 

with 



p z p z 



h p (p) = 1 - (P-j ' " . (3.20) 

Parametrising the position of the event horizon by p h , the background temperature takes the 
form 

j" = 1 ~ p C3 21) 

2tt(5-pK" V ' ' 

3.1 Generalised Conformal Symmetry 

The "near-horizon" geometry generated by a stack of N Dp-branes discussed in the previous 
section is dual to an SU(N) supersymmetric Yang-Mills theory in (p + l)-dimensions, with 
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coupling constant g Y M which is generally dimensionful, except for p = 3, and with a trivial 
RG flow which is just due to the dimensionfulness of the coupling constant. 

Following [16], let us consider a generalised version of the (p+ l)-dimensional Supersym- 
metric Yang-Mills Euclidean action 

S P + i = ~ j d^x^^-^TriF^Fn + l^^^D 2 -2Z±R% + » X \ + 

+ ^Tr{[X,X] 2 }}, 

(3.22) 

where X 1 are scalar fields, with / = 1, . . . , 9 — p, gffl is the background metric and <E> (0 ) 
is a background scalar field. The usual Supersymmetric Yang-Mills action is recovered for 
9^1 = fifiv an d $( ) = <?ym- The action (3.22) turns out to be invariant under a Weyl 
transformation of the following form 

5(0) — > e 2 "g i0) , X — > e"^!, — > A„, $ (Q) — > e^" 3 ^). (3.23) 

Defining the stress-energy tensor Ta V and the scalar operator $ as 

^ 45. f 2 5S (p + x) dcf 1 5S( P + !) , . 

one can obtain diffeomorphism and Ward identities 

V(r^) J + (O,).,^,,) = 0, (T/), 7 + (p - 3)* (0) <O,> J = 0, (3.25) 

with (•) j indicates the expectation value with respect to a source J. For g^l = 5^ u and $ (0 ) = 
<7ym t ne fi rs t equation in (3.25) provides the stress-energy conservation, while the second 
one the tracelessness of the stress-energy tensor and therefore the restoration of conformal 
symmetry. In particular, notice from the second equation of (3.25) that conformal symmetry 
is broken because of the dimensionfulness of the coupling constant. 

In a theory of this type, the entropy S at finite temperature T has to scale as 

S = c(g 2 cS (T), N, ...)U P TP, (3.26) 

with Q p being the spacial volume, g 2 a (T) as defined in (3.9) and c (g 2 B (T), N, . . .) is a generic 
function of dimensionless parameters. 

Furthermore, the two-point functions of an operator O need to have the following form 

(O(x)O(O)) =n(f {g^x), N,...) J^J , (3.27) 

with g 2 s (x) = g 2 M Nx 3 ~ p , f (g 2 B {x), N, . . .) being a function of dimensionless parameters, 
and 1Z provides the renormalised version of its argument. More precisely, these two objects 
are the same for i / 0, while at x = they differ by infinite renormalisation. 
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3.2 Holographic impurities as probe D(8 — p)-branes 

We would like now to introduce localised degrees of freedom in the (p + l)-dimensional U(N) 
SYM theory. Before looking at general values of p < 5, let us review the p = 3 case 
following [47]. 

3.2.1 Wilson loops and Db-branes 

Let us start with considering a single probe D5-brane in AdS§ x in such a way that it 
wraps an 5 3 C S 5 and extends along the time and radial direction, so that the induced 
metric on the world- volume is AdS2 x S s . This system is the decoupling limit of a system 
with N D3-branes in flat space and a single probe D5-brane, whose degrees of freedom are 
carried by the 3 — 3 strings (the "ambient" gauge theory), the 3 — 5/5 — 3 strings which 
are localised in the codimension-3 defect identified by the intersection of the two types of 
D-branes, and finally the non-dynamical 5 — 5 strings. The action for the defect (point-like) 
degrees of freedom can be obtained by performed T-duality on the D0/D8 system, obtaining 



S; 



1111 |.) 



J dt i¥d t V + (A t + Xjv 1 + at) * - na t , (3.28) 



where At and X I are respectively the time-component of the gauge field and the scalars in 
M = 4 SYM, v 1 is a unit vector, ^ is the resulting fermionic impurity field, at is the non- 
dynamical gauge field on the D5-brane, and n is the unit of background gauge field localised 
on the point-like defect. 

Integrating out the degrees of freedom associated to the D5-branes in the decoupling 
limit, it has been shown [47] that the introduction of the D5-branes corresponds to a 1/2- 
BPS Wilson loop operator in M = 4 SYM in the anti-symmetric representation of the gauge 
group U(N). 



3.2.2 General case 

As a natural generalisation of the D-brane construction reviewed above 3 , one may think to 
introduce extra degrees of freedom in the boundary theory by considering in the bulk probe 
D(8 — p)-branes wrapping an S 7 ~ p C S 8 ~ p and extending along the radial direction. The 
embedding of such branes in the Dp-brane background can be described by two functions 
x p = z(p) and 9 = 9(p), where 9 is the angular coordinate introduced in (3.16). For 
the time being, let us keep both the two embedding functions. The induced metric on the 

3 A comment is now in order. While indeed the construction we are going to consider introduces point-like 
degrees of freedom, a neat correspondence between the D(8 — p)-branes and Wilson-loop operators has not 
been shown yet. 
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world- volume of the D(8 — p)-brane is 



dsl_ p = gapdCdC? = 

= ~ hpi ^ dt2 + [1 + \{p)(z'? + uIp\(p)(0') 2 ] + 4 sin 2 9 dtf 7 _ P , 

(3.29) 

which is substantially the AdS2 x S 7 ~ p geometry with the two functions z(p) and 9(p) control- 
ling, respectively, the embedding of branes in the (p + 2)-dimensional (conformal)-j4(i5 man- 
ifold, where they wrap a (conformally^AdS^ subspace, and the embedding in the transverse 
space. The action for the probe branes is the sum of a Dirac-Born-Infeld and a Wess-Zumino 
term 

S D(8 _ p) = -T D(8 _ p) J d 9 ~ p ( e-tyj-det {g a/3 + {2^a')F a p} + T D(8 _ p) (W) j F A C 7 . p , 

(3.30) 

where g Q( g is the pullback of the near-horizon string frame metric and F is the world- volume 
abelian gauge field strength. The Wess-Zumino term in the action (3.30) acts as a source of 
the electric component Ftp of the world-volume gauge field and, therefore, one cannot put 
consistently Ftp = 0. Let us take C, a = (t, p, 6 l ) as world-volume coordinates, where the l 
(i = 1, ■ ■ ■ , 7 — p) parametrise the S 7 ~ p sphere in (3.16). Assuming that the electric field is 
independent of the angles 6 l , one can write the action as: 



Sol 



(S-p) 1 



9-p 
2 



NB P 2 n 7 - p u 7 p - p I dtdp \ [Ne^) - v p- z sm ( - p 0x 



2 
7-p 



1 + h p (z') 2 + ulp\( 



a /\2 



(2W) 2 P 4 i? p 



v (N e^>) 7 ~p _ 



1/2 



2vra' 

Br, 



(3.31) 



FtpC p {6) > , 



where 0,-j-p = 2 ir^2 E /rf^ 2 ) is the volume of the S 7 ~ p . One notices that 5 D(8 _ p) depends on 
the embedding function z and on the world- volume gauge field just through their derivatives. 
Therefore, there are two first integrals of motion 

3£ D(8 _ p) _ 3£ D(8 _ p) _ 



dFt 



const., 



(3.32) 



tp 



where £ D(8 _ p) is the Lagrangian density of the D(8 — p)-brane, whose expression can be read 
from (3.31). The first constant can be fixed through the quantisation condition [20] 



d 



7-pn d-£p(8 _ p ) 



S 7 ~p 



dF 



n 



(2W) T F1 , n G Z. 



(3.33) 
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The integer n in (3.33) represents the number of fundamental strings (quarks) dissolved in 
the probe brane. Computing the left-hand-side of (3.33), the first integral of motion related 
to the world-volume gauge-field becomes 

n sin 7 ~ p 9 Ftp 



C P (0), 



where we have introduced the rescaled world- volume gauge field F a p, defined as 

dc f 2ira' p 2 F aj3 



D 



v (Ne<t>) 7 -p 



(3.34) 



(3.35) 



and the constant Cf is defined as 



9-p 



^~D(8 - p) Bp 2 ^Ij—pUp 



7-p 



2v^F 



(2? 



(3.36) 



The equation of motion (3.34) can be easily solved for the gauge field-strength Ftp to get 



Ftp — C(p,n)(^) 



y\2 



1 + hp(z') 2 + u 2 p 2 h p (9 



sin 



2(7-p) 



1/2 



where the function Ci pn \{9) has been defined as 



C M (0) = Cp(6) + c f - 



(3.37) 



(3.38) 



As far as the first integral of motion related to the embedding function z(p) is concerned, it 
acquires the form 

p- 2 sm 7 -P 9 (N e^)~ h p z' 



[l+h p (zi) 2 + u 2 p 2 hp(6i)i - Ftp} 112 



(3.39) 



where c z is a new constant related to c z by a rescaling 



2v^ 



8-p 



9-p 



NF 



^d(8 - P )N Bp 2 f^7_ p Up p 



Fl y 



7-p\ B 



(3.40) 



Let us now evaluate the energy of the system. By performing a Legendre transform, the 
Hamiltonian H of the D(8 — p)-brane is given by 



H 



57-p 



dp 



Ft 



tp 



'D(8 - p) 



dF t 



tp 



(8-p) 



(3.41) 



By using the explicit expression of the Lagrange density written in (3.31), one gets 

dp (N e+) ^ p~ 2 [1 + hp(z') 2 + p%(9') 2 ] 1/2 x 

1/2 



H = F u{8 _ p) NBp 2P n 7 ^ p u 7 -f 



(3.42) 



sin 



2{7 - p) + C 2 M (9) 
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Let us now discuss some possible configurations for the D(8 — p)-brane. In particular, we 
will focus on the case in which the position of the probe branes in both the (conformally) 
AdS space and in the transverse space is fixed (straight flux tubes with 9' = z' = 0), and 
on the configurations in which the probe branes are allowed to bend either in the extended 
directions or in the transverse space. 

4 Straight Flux Tube Configurations 

Let us start with the analysis of the configurations in which both the angular and linear co- 
ordinates 6 and z are constants. Notice that these configurations introduce a localised defect 
on the boundary theory. Following [20], the stable configurations are the ones minimising 
the energy (3. 42), i.e. the ones that satisfy the condition: 

2 



dH 



= const 
z — const 



T D{8 _ p) NByn 7 _ p u 7 ~P I dp (Ne*) 7 ~ P p~ 2 x 



x (7-p) sin 7 p i 
where the function A p n (#) is defined as 



Ap,»(0) 



(4.1) 



Sin 2(7- P ) Q + C 2 (0)1 V2 ' 



Ap, n (0) = sin 6 -^cosfl-C p , n (fl). (4.2) 

The least energy condition (4.1) can be satisfied if and only if 

sine = 0, or A p ,„(6») = 0. (4.3) 

In the first case, such configurations occur at 6 = 0, ir which are points in which the sphere 
S 7 ~' p shrinks to zero size. In the second case, instead, the branes are located at 6 = 0( p ^ n ) 
with #( Pjn ) defined by the condition itself: A Pjn (0(p in \) = 0. The functions A p ^ n (9) for different 
values of p are listed in appendix A. As it is clear from (3.38) they depend on the quantisation 
integer n and on the rank of the gauge group through the combination 

v = — , (4.4) 

which we will refer to as the filling fraction. The reason for this name is the fact that, for a 
given value of p, only for 1 < n < N there exists a unique solution for the angles 0(p,n) m the 
range < 9(v,n) < 7r - The integer n represents the number of impurity fermions introduced in 
the boundary theory and the ratio v takes values in the range < v < 1. The upper bound 
of n is a manifestation of the so-called stringy exclusion principle and is a piece of evidence 
supporting the identification of these brane configurations as the holographic duals of Wilson 
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lines in the antisymmetric representation of the gauge group. Indeed, in the conformal case 
p = 3 this identification was explicitly checked in refs. [47,48]. 

Prom the explicit formulas of the K p ^ n {9) functions displayed in appendix A, one can 
verify that the angles 9( p , n ) satisfy the relation: 

6{p,N-n) = 7T - 6>( pjn ) . (4.5) 

It follows that changing the polar angle of the embedding as 9 — > tt — 9 is equivalent to the 
particle- hole transformation v — > 1 — v. The energy density (tension) of these configurations 
was derived in [20] and it turns out to be 

NT r ( 7 -f) 



2 

It follows from (4.5) and (4.6) that the tension satisfies the relation: 

£(p,n) £(p,N —n) i 

(4.7) 

i.e. it is invariant under the transformation v — >■ 1 — v. 

The electric world-volume field Ftp for the flux-tube configuration 9 = 0(p^ n \ can be 
obtained from (3.37) by using that C( p ^(9( p ^) = (sin 6 T ( pn )) 6_p cos#( Pjn ). One gets: 

T t p = cos ^p^) . (4.8) 

Notice that Ttp changes its sign under the transformation v — > 1 — v. 

As is was proven in [20], in the case of the zero temperature background, the flux-tube 
configurations described above are the solution of a first-order BPS equation for the embed- 
ding and the world-volume gauge field of the probe brane. One can show that this equation 
is the one that is obtained by imposing kappa symmetry to the D(8 — p)-brane in such a way 
that it preserves 1/4 of the super symmetry. 

4.1 Impurity Entropy 

Let us investigate some aspects of the thermodynamics of this class of systems. First of all, 
for p < 5 the relation (3.21) between the temperature T and the position of the horizon p 
can be inverted to obtain 

- ^w=w- (4 - 9) 

In order to compute the free energy and entropy of the flux tube configuration, let us evaluate 
the Euclidean action of one of such configurations that extends from the horizon at p = p h 
until a cutoff value of the radial coordinate p = e. If we define f3 p as 
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we get 



D(8-p) 



f -2- 

_ °(P>") o n 5-p 

— rp Hp C> V Ph 
on-shcll ± 



Ph 7-p 

dp p 5-p 



(4.11) 



To arrive at the expression (4.11) we have integrated over a periodic Euclidean time circle of 
period 1/T. Notice that the integral (4.11) diverges as e — > and therefore it needs to get 
renormalised in a neighbourhood of the boundary. It easy to see that it is renormalised just 
by a term proportional to the volume of the boundary, regularised with p = e 

5 — p „) 



7-Eucl 
'D(8-p) 



B p [Ne* 



2 

7-p 



Ett\ e 



(4.12) 



where g« | £ = e 2 is the tt component of the induced metric (3.29) at p = e. Therefore, the 
renormalised action is given by 



Eucl 
D(8-p) 



5-p S t 



(P, n) 
rp Hp 



ftp B p 



(4.13) 



Let us rewrite (4.13) in a more convenient form. By using the definitions of j3 p and B p written 
in (4.10) and (3.11), one easily proves that: 



(3 p Bp 



2v^F 



(5-p)T P 



Eg) 

(7-p) 2 IT 



1 

5-p 



(4.14) 



where g cS (T) is the gauge theory effective coupling constant at the temperature T, defined 
as 

£ S {T) = gl M NT^ 3 . (4.15) 

By using the explicit expression of £( Pi „) written in (4.6) one can recast the renormalised 
action (4.13) as 

(4.16) 



rBucl 
J D(8-p) 



,iV(sin^ n) ) 6 p [g 2 ce (T)Y'-v 



where the coefficient c p is given by: 



(i — p 
5—p 



(4.17) 



From the renormalised action (4.16) it is immediate to extract the free energy for the straight 
flux tubes, namely 

= -c p N ( sm6 M ) 6 " p [g 2 a (T)] ^ T . (4.18) 



™tr 

r D(8 - p) 



rjs tEucI 

1 1 D(8-P) 



If we ignore non-abelian interactions among the D(8 — p)-branes, the impurity entropy gen- 
erated by M D(8 — p)-branes is just M times the result obtained for M = 1, namely 



°D(8-p) 



-M- 



dR 



D(8 - p) 



dT 



1c p 
5 — p 



MN {sm6 M ) 6 p [g 2 a (T)}^ 



(4.19) 
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It is interesting to study what happens in case of small filling-fractions v = n/N. The impu- 
rity entropy (4.19) depends on the filling-fraction v through the factor sin 6-p (9 (p n) contained 
in the energy density £( p , n )- In the small filling-fraction limit the angle Q(p t n) ls a ^ so small 
and one can check [20] that, at leading order in n/N, one has: 



6— p 



2^ 



n 



r{ 7^,\ N 



H 



(4.20) 



Then, the expansion of £ (p , n) is given by 



lim £, 



— ->o 



(/>. 



nT F1 



+ 



6 — p 



P 2(8 -p) 

and, as a consequence, the impurity entropy can be written as 



2^F 



r 



(5? 



2 

6— p 



(4.21) 



lim 5 



D(8-p) 



r? 



with a 



Ma p [ 5e 2 ff (T)] 

def 4^/? 



1 



W \ 6-p 



+ 



P 



(7- P ) 2 7T 



(4.22) 



One can verify that (4.19) and (4.22) for p = 3 coincide with the impurity entropy computed 
in [23, 24] for the holographic dual of the maximally supersymmetric Kondo model. Our 
results generalise those in [23, 24] for non-conformal Dp-brane backgrounds. 

Prom (4.22) we notice that the dependence of S on the filling fraction is, in general, 
non-analytic. However, in the particular case of the D4/D4 system, the impurity entropy 
turns out to be analytic in the filling fraction. Actually, it is possible to obtain a simple 
closed expression of S£l as a function of v. In particular, for such a case the system is in the 
minimal-energy configuration if the probe D4-brane is located at 9 = # (4i n) such that 



cos 9 



(4,n) 



1 - 2za 



(4.23) 



satisfying the condition (4.3). Using (4.23) the expression for the entropy density for p = 4 
can be written as 



2nM 
9 



9 2 AT) (1 - v) 



(4.24) 



where g^ s (T) = g\ u NT in this p = 4 case. 

The internal energy £ is obtained from the free energy F by means of the standard 
thermodynamic formula 

£ = F + TS . (4.25) 
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For the straight brane configuration one can immediately obtain £ by combining (4.18) and 
(4.19), namely 



£d(8- p) = H^^D(8- P) • (4-26) 



p — 3 

T 

Notice that £^ vanishes, while £|v r 8 ^ becomes negative for p < 3. 

In order to complete our thermodynamics analysis, we can compute the impurity specific 
heat from the impurity entropy, namely (4.19) 

C&-„ = = ^Z|a MiV (sin^ n) )^ [fli(r)] ^ , (4.27) 

which vanishes for p = 3 and is negative for p < 3, giving a signature of a thermodynamic 
instability. From (4.27), it is easy to compute the impurity susceptibility 

i 



dT 

which is constant for p = 4 



XnV.) « — ^ <* (P-3) 2 L ^7 J = (P-3) 2 L ^^ , (4-28) 



5 Hanging Flux Tube Configurations 

Let us now consider the configuration in which 9 is a constant and the embedding of the 
probe brane in the (p + 2)-dimensional (conformally)-yl(i5 manifold is described through the 
scalar z(p). The probe branes are located at 9 = 9( p , n )- 

Analysing the form of the first integral of motion (3.39), one can notice that, in the case 
of black hole embedding phase, the regularity condition at the horizon fixes the (rescaled) 
first integral of motion c z to be zero, and, as a consequence, the embedding function z must 
have a trivial profile. 

Focusing instead on the configuration in which the probe brane lies completely outside 
the black hole, (3.39) implies the existence of a turning point. Inverting (3.39) with respect 
to z' and integrating the obtained equation, one gets 

<p) = ± / = 9 c ( r ] = . Z z p - n • (5.i) 

sin v 



p 



^ p [(Ne^{p<)-^ p -(c { i' n) ) 



>p,n 



It is easy to see from (5.1) and (3.40) that the constant ci P ' n) is actually related to the energy 
density (4.6) 

c ( P , n) = —^—^ (5.2) 

c (p,n) a V 

where c z is the constant introduced in (3.32) and the position p t of the turning point is given 
by 

Pt = — -5^~ » ( 5 - 3 ) 
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Let us next introduce a new radial variable a, related to p by means of the expression 

P 



a 



Ph 



Clearly, the turning point in this new variable is just 



2(7-p) 



(5.4) 



(5.5) 



Moreover, the separation on the boundary between the re-quark and the re-antiquark is given 
by 

o 2 (P. n) 7-p 

L = at P X{a t ) , (5.6) 



where I(a t ) is the integral 



da 



7-P 

a 5 ~p 



7-p 



,7-p 



>Z^£ 



1/2 • 



(5.7) 



By using (5.5) one can eliminate in (5.6) the constant ci P in favour of a t , with the result 



L = 2 Ph A/1- (aO^ X(ot) 



(5.8) 



By applying the same techniques as in appendix B, the integral (5.7) can be obtained in 
terms of the hypergeometric function. One gets: 



2fo) 



1 7-p 



5— p 



VTCT t F( 1 1 



1 3 1 2|=£ 



2 1 7-p ' 2 7-p 



; <7 t 



(5.9) 



,2(7-p) / 

Let us use these results to write L in a closed form. Actually, it is very convenient to define 
the parameter 7 as 

7 = af 1 ^ . (5.10) 
It follows from this definition that < 7 < 1. Moreover, L can be written as 



L 



7 — p v 7 ~p y „, 5 ,r p 



(5-p)^rT r 



2(7=^ 



p ' 2 7 — p 



where we used (3.21) to eliminate ph and we wrote the result in terms of the temperature T. 
The energy for this hanging configuration is 



E — 2£(p ,«) PpBpCTt p 



da V 1 - <t 2 ^ 



7-p 



7-P 



(5.12) 
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where the cut-off e has been introduced to regularise the Hamiltonian near the boundary 
<7 = 0, where the integral diverges. It can be regulated by minimal subtraction, i.e by 
the energy corresponding to a disconnected configuration. Therefore, let us look at the 
divergences near the boundary a = 0, where the embedding function z(p) and its first 
derivative have the following asymptotic expansion 



z'(a) 
z{a) 



1=2. 1 , 

a s- P 1+ [ 1 + CTt 



L 5 — v pi , s <2<>=p 



Pi 



B-f \ a 2 5=t + (o- 4 5. 



T 2 ± 2(6 -p) P P C 



cr z 5-p +0 (a 2 5 ~p 



13 -2p 



(5.13) 



It is easy to read off the divergent boundary Hamiltonian 



If 



(5 - p)£ {Pin) (3 p B p e 5 -p. 



(5.14) 



As a consequence, one needs to add just a counterterm which is proportional to the volume 
of the boundary, regularised with a = e 



H 



-(5-p)£(p in) Bp [Ne 



2 

7-p 



Stt\ 



(5.15) 



As noticed in [27], the divergent Hamiltonian (5.14) is exactly the divergent boundary Hamil- 
tonian that one would obtain from probe D-branes in an ^4<lS-space of dimensions 2/(5— p)+l, 
and the counterterm (5.15) is actually given in terms of the volume of the boundary of this 
AdS-sp&ce. 

Therefore, the renormalised Hamiltonian is given by 



H 



and the free energy can be written as 



r T>(6-p) 



2 £(p. n) ftp Bp 



da 

7-P 
(J5-P 



lim 



a 



H 



+ H 



(5.16) 



7-P 
5— p 



2 7-P 
l — O" 5 ~P 



2 7-p 

a t 5 ~ p - a~-v 



1 



>2=p. 



(5.17) 



Let us now change to a new variable £ in the integral (5.17), defined as 

e = 



a 



2 

5— p 



(5.18) 



After this change of variable, it is easy to verify that F^, 8 _ > becomes 



r D(8-p) 



7-P 



d( 



7-p 



VC~ P - 1 



(5.19) 
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where 7 is related to a t as in (5.10). The term inside the square brackets in (5.19) is just the 
integral J(7 — p, 7) calculated in appendix B. Thus, we can write 



r D(8-p) 



(5-p)£(p,„) PpBp-f r-p J(7_p 57 ) . 



(5.20) 



The prefactor in (5.20) can be related to the free energy of the straight brane. Indeed, 
according to (4.13) and (4.18) one has (5 — p)£ {Pin) (3 p B p = F^L_ p y Moreover, by using the 
explicit value of the integral J(7 — p, 7) (see eq. (B.6)), one gets 



r D(S-p) 



9 F BtI <"V 7- 



r( ^2 

1 7-P 



r 



5— p 
2(7-p) 



1 



5 — p 



2 ' 7-p ' 2(7 -p) 



■ 7 



(5.21) 



The constant 7 appearing on the expressions of the length L and free energy -Fp (8 _ p) (eqs. 
(5.11) and (5.21)) parametrises the turning point of the hanging brane configuration, i.e. 
how much the brane penetrates into the bulk. The configuration with 7—7-0 corresponds 
to p t — > 0, which means that the brane lies completely on the boundary in this limit. On 
the contrary, when 7=1 one has p t = p\ t and, thus, this configuration corresponds to two 
disconnected straight branes ending on the black hole horizon. It is easy to check that our 
formula (5.21) is consistent with this interpretation. Indeed, as 



1 



5 — p 



2 ' 7-p' 2(7- p) ' 



1 \2(7-p) 



7T W tj— p 
1 7-p 



(5.22) 



one gets 



* r D(8 -p) ■ 



(5.23) 



Finally, we can now compute the one-point correlator of the boundary operator O z as- 
sociated with the scalar function z. In order to perform this computation, we consider the 
Euclidean version of the action (3.31) at 6 = 9( P:n ) 



7-Eucl 
1 D(8-p) 



■D(8-p) 



dtdp [Ne 



7-P ^2 
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' (p, n) ■ 



1 + h p (z') 



sin 6 p 9( Pi „) 



•1 1 



tp 



1/2 



2(2vr)P- 2 
1/2 



P-3 

5_P I d p \ (5-p)(7-p) 
Un 



(5.24) 



dtdp p 5 ~p<sin p 9( p n) x 



i + h p (z'Y-^ p +F tp c p (e ip , n) ) 



where the second equality has been written using the explicit expression for the dilaton 
(3.15) and U p (3.7) as well as the relation (3.36) and the expression for the energy density 
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(4.6). It is interesting to notice that the overall term p~~^-p in (5.24) can be seen as the 
determinant of an AdS black- hole metric from two-dimension to 1 + q, with q = 2/(5 — p). 
Therefore, we can think that our probe D(8 — p)-brane configuration is actually in a higher 
dimensional AdS black-hole background (times S s ~ p ) and that the induced metric on the 
brane has an (asymptotically) AdS\ +q factor, with q integer - then we can perform the 
analytic continuation to the value 2/(5 — p). This action is renormalised exactly by the same 
term (5.15) used to renormalised the Hamiltonian - notice that if one writes the explicit 
expression for the dilaton in (5.15), the counterterm can be seen as proportional to the 
volume of the boundary of AdS\+q. In this way, we can compute the one-point function for 
the operator O z dual to the embedding function z(p) as 



lim 

e->0 e 



r 7-Euci 

J D(8-p)| 



A, 



9dAdS 1+q | 



5z(e) 



(5.25) 



where gaAds 1+q \ € is the determinant of the metric at the boundary of AdS 1 + q regularized at 
p = e, and A z = 2/(5 —p). Using equation (5.25), the explicit expression for the one-point 
correlator is 



n) Bp 



sin 6 p 8 (p , n) 
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; sin# 
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e5_p \ g 9 Ad Sl 



+i\e 



i + % {z'Y 



1/2 



(5.26) 



where the last equality has been obtained by using the equation of motion (3.39) for z 4 , with 
c z has been defined in eq (3.32). 



6 Dimer Thermodynamics 

In this section we will study in detail the thermodynamic properties of the connected config- 
uration of the D(8 — p)-brane. First of all, we will rewrite the free energy (5.21) in a more 
convenient way. With this purpose, let us define two functions h(-j) and g(-j) as follows: 

, , . def - 1 „/ 1 11 1 \ 

h(j) = 7 r- P F{ - - , ; - - ; 7 , 

V 2 7 — pi 7 — p J 

, . d ef I L. / „/ 1 13 1 \ , 

g(i) = V 7 - VT^f(-,i- — ; _-_ ;7 j , (6.1) 

4 One has to remember that this computation has been performed by using the Euclidean action, as usual. 
If one repeats the same calculation using the Minkowski signature, the result would be the same up to a sign. 
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as well as the constant A„, given by 



-v 

A P = y~ P \ • (6.2) 

5— p 



Then, the free energy (5.21) can be written as 

F^s- P) (T) = 2A P F*%_ p) (T) h( 7 ) , (6.3) 

where we have stressed the fact that both F^ g _ p) and Fg%_ * depend on T. The dependence 
of F^l_ p) on T is shown explicitly in (4.18). Moreover, F^ s _ p) depends implicitly on T 
through its dependence on the parameter 7. Actually, by rewriting (5.11) as 

A p5 (7) , (6.4) 



(5 — p) ir L 

we get the explicit relation between the temperature T and the parameter 7 (for fixed length 
L). In order to study the dependence of the temperature on 7, let us define the dimensionless 
reduced temperature T as 

r = |4ir. (6.5) 

(7-p)A p 

It follows from (6.4) that T(^f) is just given by the function 5(7), namely 

T( 7 ) = 9(l) • (6.6) 

Let us now analyse the behaviour of 7~(7) near 7 ~ 0, 1. From the explicit expression of 5(7) 
in terms of the hypergeometric function (see the second equation in (6.1)), one can verify 
that 

Tin) ~ T^" 7 ^ , 7~ , 



r(7) « K \ 7 "\ V^l log(l - 7) , 7 « 1 , (6-7) 

^r(i-^) 

and, therefore, 7"(7) — > as 7 — > 0, 1. Actually, by plotting the function T{^) for different 
values of p one discovers that it reaches a maximum at some value 7 = 7™ ax with < 7 p nax < 1. 
These plots are shown in figure 2 for < p < 4. 

The fact that T has an upper bound implies that above a certain temperature 7^, max = 
7"(7 = 7 p nax ) (which depends on p), only the disconnected solution exists. The values of 
7™ ax and 7^, max can be obtained by requiring the vanishing of the first derivative of 5(7) • We 
can calculate g' (7) by differentiating the right-hand-side of the second equation in (6.1). By 
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Figure 2: On the left we plot the reduced temperature T versus the parameter 7 for < p < 4. 
As p is increased the maximum of the curve is shifted towards lower values of 7. On the right 
we plot Af p versus 7 for < p < 4. The points where the different curves cut the horizontal 
axis define the 7 parameter 7* of the dimerisation transition . 



doing so we get a term which contains the derivative on the hypergeometric function, which 
we simplify by using the relation 

7(1-7) A. F(a,6;c; 7 ) = (c - a) F(a - 1, b; c; 7) + (0 - c + 67) F(a, 6; c; 7) . (6.8) 

0(7 

Then, g'ip/) can be written as a combination of two hypergeometric functions which, remark- 
ably, can be simplified to be: 

Therefore, the maximum of the reduced temperature T(j) occurs at 7 = 7^ ax , where 7™ ax is 
the solution of the equation: 

F {-i l -rh^-rh ii! ")- - (6 ' 10) 

The numerical values of 7™ ax for < p < 4 and of the corresponding reduced temperatures 
-y-max j^yg Deen written in table 1. One notices that 7™ ax becomes smaller as p is increased, 
while 7^ max grows with p. 

For T < 71 max it is clear from figure 2 that there are two connected solutions with different 
7 for a given value of T ■ Recall that 7 parametrizes the holographic coordinate of the turning 
point of the connected configurations. In addition, we have the disconnected configuration, 
which is the only one that exists for T > X, max and is competing with the two connected 
ones when the temperature is lowered below 7^ max . In order to determine which one of these 
three configurations is thermodynamically more stable we have to find out which one has the 
lowest free energy. One can establish numerically that the configuration with higher 7 (i.e. 
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p 


max 

lp 


srmax 
lp 


* 

7 P 


lp 


Ae* 





0.6141 


0.6916 


0.2745 


0.5927 


1.5150 


1 


0.5965 


0.6979 


0.2557 


0.5995 


1.5274 


2 


0.5693 


0.7078 


0.2285 


0.6101 


1.5465 


3 


0.5215 


0.7254 


0.1855 


0.6293 


1.5799 


4 


0.4155 


0.7658 


0.1102 


0.6752 


1.6541 



Table 1: Numerical values of 7p iax , 7^ max , j p and T p * for < p < 4. On the last column we 
have displayed the reduced latent heat of the dimerisation transition, as defined in (6.38). 



extending deeper in the bulk) has higher free energy than the connected one with lower 7. 
To compare the latter with the disconnected one, let us define the reduced difference of free 
energies Af p as: 



def 1 D(8 - p) 



9 pstr 



9 Z?str 
Z f D(8-p) 



(6.11) 



From (6.3) we get that, as a function of 7, Af p is given by: 

Af p = 1 - A p h( 7 ) . 
Clearly, by its definition (6.11) and by (5.23) it follows that Af p 

1 

near 7 ~ one can check that Af p diverges as Af p ~ — Ay, 7 



7-p 



(6.12) 

as 7 — > 1. Moreover, 
-00. The functions Af p 
have been plotted in figure 2. It turns out that Af p (-j) vanishes for an intermediate value 
7 = 7* which, according to (6.12), satisfies: 

Kil) = 4-, il < 1 • (6-13) 



The different numerical values of 7* for < p < 4 have been collected in table 1 , together 
with the corresponding values of the temperature T* = T(7*). As shown in table 1, 7* < 
7™ ax , which means that the connected configuration with lower 7 is the more stable one 
when T < T p * ■ At T = T p it takes place a dimerisation transition of the type advocated 
in [21], in which two disconnected branes are recombined to form a hanging brane and the 
impurities of the boundary theory condense to form bonds or dimers. In order to complete 
the thermodynamic description of this dimerised phase, let us calculate the entropy of the 
connected configuration. By computing the derivative of (6.3) with respect to T, we obtain: 







dT 

By using that (see (4.18) and (4.19)) 



S^s- P) (T) = -T^ F- 8 _ p) (T) = 2A P 5* _ p) %) - F*- r) ^ {^) ■ (6-14) 



dh (dT_y l 



= ~^TSSl s _ p) . ( 6 - 15 ) 
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together with (6.4), the relation (6.14) can be written as 

)(T) = 2A P 5n ( r 8 _ p) 



°D(8 - p)\ 



M7J + — s~ 5(7) -7: 



(6.16) 



2 ^" 5 '( 7 ) 

The derivative of g(i) has been calculated above (see (6.9)). We now compute h'(i). To 
carry out this calculation we will differentiate the right-hand-side of the first equation in (6.1) 
and we will use the relation: 



F(a,b;c;i) = — F(a + 1, b + 1; c + 1; 7) . 

cry c 



(6.17) 



Proceeding in this way, one can represent h'(i) as a combination of two hyper geometric 
functions which, remarkably, can be written as a single hypergeometric function. One gets 



1 



7 — p 



1 3_ 

1 7-p _p 



1 



1 



1 1 



7 — p 2 7 — p 



■ 1 



(6.18) 



Amazingly, the hypergeometric functions on the right-hand-sides of (6.18) and (6.9) are the 
same and, therefore, the ratio h'{i)/g'{i) is remarkably simple, namely 



1-7 



g' (7) 5-p 1 7 

Therefore, the entropy of the connected configuration can be written as 



'D(S-p) 



(T) = 2A P SSI 



8-p) 



Hi) 



1-7 



5(7) 



(6.19) 



(6.20) 



Again, the magic of the hypergeometric functions allows us to write the term in brackets in 
(6.20) in a simplified form. One can actually check that 



Hi) 



1-7 2 (6-p)(7-p) , 

9(V = 777. iTTTi \ a (V > 
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(19-3p)(5-p) 



.21) 



where 17(7) is a new function defined as 

ad)^! 1 '^ ^(1,1-^1 



1 



7 



p 2 7 — p 

By using this result, one can easily show that the entropy can be written as 

S^ s _ p) (T) = A p S** a _ p) (T)a(i) , 
where the coefficient A p is given by: 



(6.22) 



(6.23) 



A., 



r 2- * 

def V 7 ~P 



r 



7-p 



7T . 



(6.24) 
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As a check of the equation (6.23), let us take the limit 7 — > 1. One can easily show that 

2 



(7(7 = 1) = j- , 



and, therefore, one has 



lim S D V p) (T) 



2^V P) (T) 



(6.25) 



(6.26) 



as expected. In order to study in more detail the entropy of the dimers, let us define the 
reduced entropy s p (T) as the ratio 



^D(8 -p) 



(T) 



2 CO 



(6.27) 



It follows from (6.23) that s p (T) is given by 



A, 



(7(7) . 
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By plotting the right-hand-side of (6.28) versus 7 one realises that s p is a monotonic function 
of 7 which grows from zero to one as 7 varies in the interval [0, 1]. Let us study the behaviour 
of the entropy of the dimer near T ~ 0, which corresponds to 7 ~ 0. From the definition of 
(7(7) in (6.22) it follows that £7(7) ~ 7 7 —p — > for small 7. Moreover, it is straightforward 

2(7-p) 

to prove from the expression of the reduced temperature in (6.6) that 7 ~ (7J 5 ~ p for 
small 7. Thus, we have for small T the reduced entropy of the dimer is given by 



-f (r) 5 - 



(6.29) 



Let us now use this result to obtain the low temperature behaviour of the entropy S- D ^ g _ p y 
By using (4.19) for M = 1 and (6.5), it follows that, for low T, the entropy of the dimer of 
length L behaves as 
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where a„ is a numerical coefficient given by: 
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Notice that we arranged the right-hand-side of (6.30) in such a way that all factors are di- 
mensionless. The first two factors in parenthesis contain the dependence on the filling faction 
and on the Yang-Mills coupling respectively, while the last factor contains the dependence on 
the temperature. Notice that the power of T in (6.30) is positive for all values of p (contrary 
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to what happens with the straight-brane configuration) which means that the entropy of the 
dimer vanishes at T = 0. In particular, for p = 3, the entropy of the dimer configuration of 
D5-branes in the AdS*, x S 5 background is given by 

N( sin£ (3in) ) 3 V\(LTf , 



'US 



(6.32) 



where the relation between the angle #(3 >n ) and the filling fraction has been written in (A. 4), 
A is the 't Hooft coupling of the bulk J\f = 4 theory and the numerical coefficient is just 03. 
Let us next compute the internal energy for the connected configuration, which is given 

by: 

^D(S-p) = Fd(8-p) + ^>D(8-p) • (6.33) 

By combining (6.3), (6.15) and (6.20), one gets 
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Let us rewrite (6.34) by using (6.21) to eliminate h. After some simplifications one has 
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Notice that, due to (6.25), the right-hand-side of (6.35) reduces to 2£p ( r 8 _ p) (T) as 7 — >■ 1, 
as it should. In order to characterise the dimerisation transition, let us define the reduced 



latent heat Ae p as 



Taking into account that 2£^J 
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(6.36) 



3.37) 



By construction Ae p (7) — > 1 when 7—7-1. Moreover, Ae p (7) — > +00 as 7 — > since 

1 

Ae p (7) ~ (5 — p) A p 7 7 -p for small 7. By plotting Ae p (7) as a function of 7 for different 
values of p one can verify that Ae p (7) > for 7 £ [0, 1]. It follows that the reduced latent 
heat at the dimerisation transition, given by 



Ae; = Ae p ( 7 = 7 ;) 



(6.38) 



is always positive and therefore the dimerisation transition is a first-order phase transition. 
The values of Ae* for different p's have been written in the last column of table 1. 

We now compute the specific heat for the dimer configuration. Proceeding as in the 
calculation of the entropy, we arrive at 



°D(8-P) 



(T) = T-^S^ 8 „ p) (T) 
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(6.39) 
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Figure 3: Reduced specific heats 0^(7) for < p < 4. Only a portion of the region < 7 < 
7™ ax , where c p > 0, has been represented in the plot. The curve that grows faster (slower) 
corresponds to p = 4 (p = 0). 

By using identities satisfied by the hypergeometric functions, the derivative of 17(7) with 
respect to 7 can be simply written as 

y (7) = |z£ 7 -AiP(i, 2 - 7 i- i |- 7 i- i7 ). (0.40) 

Moreover, since T can be written in terms of the function g (see (6.4)), equation (6.39) can 
be recast as 

C D U (8 _ p) (D = A p SZ _ p) (T) [ a( 7 ) + 5 ( 7 ) ^ ] , (6.41) 

where 5' (7) and <r'(7) are written in (6.9) and (6.40) respectively. In figure 3 we plot the 
reduced specific heat, defined as 

- - C ° V " - P - 3 -(7) + 9 (7)^, (0-42) 



A P SJ%_ P) 5-p g'h) 

as a function of 7 for different values of p. Notice that c p — > 00 as 7 — > 7^ ax , since the 
derivative of 5(7) vanishes at this point. Actually, one can check that 0^(7) is positive in the 
region < 7 < 7^ iax for all values of p < 4. Notice that the denominator A p S^L _ p) in (6.42) 
is always positive and, thus, Cn (g _ p) > for 7 € [0,7™ ax ). The behaviour of C^ 8 _ p) (T) for 

9-p 

low T can be obtained directly from (6.30). It follows that Cp (g _ p) (T) ~ T 5 ~p — > as T — > 0. 
7 Impurity Fluctuations 

In this section, following [20,25] we discuss the fluctuations on the probe branes around the 
straight flux tube configurations in the zero temperature background. The goal is to check 
the stability of these configurations and to carry out the holographic renormalisation program 
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for the corresponding action of the fluctuations and to compute the one-point and two-point 
functions for the fluctuation operators. 

The action that governs the fluctuations of the D(8 — p)-brane probe can be obtained by 
expanding its DBI action around the configuration with z' = and 9 = 9( Pin )- The induced 
metric on the world-volume for this configuration at zero temperature is just obtained by 
taking z' = 9' = and hp = 1 in (3.29), namely 

C dC dC? = '^J; ^ + ul sin 2 9 {p , n) dtf 7 _ p . (7.1) 

Notice that this metric is of the type AdS2 x S 7 ~ p , with u 2 sin 2 #( Pi „) being the radius of the 
S p sphere. 

Recall that our probe branes have a non-trivial world- volume gauge field Tt p = cos #(p jn ) 
(see (4.8)). Due to this, the metric that is relevant for the fluctuations is not the induced 
metric g (0) but the so-called open string metric, which includes the effects of the world- volume 
gauge field. In order to define this metric let us consider the matrix (<? (0) + T\ , which is 
the inverse of the matrix appearing in the DBI Lagrangian for the unperturbed configuration. 
Let us write this inverse as the sum of a symmetric and an antisymmetric matrix 

'g (0) + TY 1 = g- 1 + J , (7.2) 



with J being the antisymmetric component. Then, the symmetric matrix Q a p is defined as 
the open string metric. In our case it is straightforward to obtain the expression of Q a p. One 
gets 



G aP (K a dC P = sin' 



(p, «) 



-dt 2 + dp 2 



+ u 2 dVL 2 _ 



p2 1 P ^" U i-P 



(7.3) 



By analysing the possible fluctuations of the D(8 — p)-brane probe one can identify two 
decoupled channels: 

• fluctuations x °f the Cartesian coordinates; 

• fluctuations (£, /) of the embedding function 9 and of the world-volume gauge field. 
These two channels will be studied separately in the next two subsections. 

7.1 Fluctuations of the Cartesian Coordinates 

Let us assume that the Cartesian coordinates are not fixed but, instead, they fluctuate, and 
let x be the corresponding fluctuation. In this fluctuating configuration, the induced metric 
on the probe D(8 — p)-branes becomes 

S 2 _ ~(0) , ~(i) j/-aj/-B ~(l) dcf 



dsj- P = 9^dCd(P + ~g^dCd(P, g% = p~ 2 (d aX ) (d pX ) , (7-4) 
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with being the AdS2 x S 7 ~ p zero-order metric (7.1). By plugging the metric (7.4) into 
the DBI action and by expanding to second order in the fluctuation, one gets that the action 
for x can be written as 



S 



(x) 



£(p, n) Bp sin 2 Q( p n ) 



n 



7-p 



-p' 2 Q af} (d aX ) (d 0X ) + O( X A ) 



(7.5) 



where gAds 2 an d 9s 7 - p are respectively the determinant of the AdS2 metric and of the unit 
sphere S 7 ~ p , while Q a p is the open string metric (7.3). In (7.5) we have written the global 
coefficient in terms of the flux-tube tension £( p ,„) and of the constant B p defined in (3.11). 

It is convenient to rewrite the dilaton by factorising the constant contained in its definition 
(3.15) and to introduce a rescaled dilaton tp as 
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e 2 ( 5-p) v _ N e <t>* e 2( 5 -p) v s (7. 6 ) 



so that the new dilaton 93 is just 99 = log p 1 and, in the second step of (7.6), we have defined 
a new constant <fi*. In terms of ip the action acquires the form 



-T, 



(p. 1) 



^p- 2 G a HdaX)(dpx) + 0{x A 



(7.7) 



where the effective tension r (Pi n) is given by: 



def B p (Ne**y-> 
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(7.8) 



By using the explicit values of B p , (ft* and £( p , n ) it is straightforward to verify that T( p?n ) can 
be written in terms of field theory quantities as 
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(7.9) 

It is interesting to notice that the action (7.7) is equivalent to the action of the Cartesian 
fluctuations for a brane wrapping an ^4dS'i+ -space rather than AdS% 



q(x) 



-?D(7+a-p) j d 8+<7 P ( ' — gAdS a + 1 \J 9s 7 - 
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P ~ 2 G ap (d aX ) {d P x) + 0{x q 



up to a compactification on a T a torus in such a way that the metric is 



dsi 



-dt 2 + dp 2 



P-3 



a i/-b 



(7.10) 
(7.11) 
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where a, b run on the a— 1 extra directions, and the the extra coordinates £ a take values in the 
interval [0, 2irR T ]. Notice that the metric (7.11) corresponds to the one of the AdSi +a space 
only when the prefactor multiplying the line element of the extra coordinates is e 2 ^ = p~ 2 . 
This happens when a is analytically continued to take the factional value a = 2/(5 — p). 
Moreover, the compactification radius R T is related to the tensions T D(7 + a _ v) and r (p . n) by 
means of the equation 

T D{7 + CT „ p) (27ri? T ) <J ~ 1 = Tfan). (7-12) 

This action can be considered as the action of a D(7 + a — p)-brane in an AdSi+ p + a x S s ~ p 
wrapping an AdS\+ a x S 7 ~ p submanifold. 

In order to easily compute the asymptotic behaviour of the fluctuating field x, one can 
consider the factor p~~ 2 in the square brackets of (7.10) as a further dilaton factor and then 
the action (7.10) can be conveniently written as 

S!ft r _ p = -T D(7 + CT _ P) | d*-v + °C,e 2 ^-g AdSl+a ^- p \q<# (d aX ) (d pX ) + O ( X 4 ) 

(7.13) 

As before, one can consider the dilaton factor as further extra-dimensions in the world- volume 



$7+ q - P — —T D (7 + q- P ) I d 8+q p ( ■sZ—gAds^Vds 7 -" 
where q is given by 



\G aP (d aX ) {dpx) + 0{ X A 



, (7.14) 



q = 2(6-p)/(5-p) , (7.15) 

and the compactification on the T q ~ p torus is determined by the relation between the tensions 
of the D(7 + q — p)-brane and T (j , n) , namely 

T D(7 + q _ p) (2irR T ) q ~ 1 = r CPin) . (7.16) 

Let us next consider an ansatz for the fluctuation field x such that is factorises as the 
product of functions of the AdS space and of the S 7 ~ p sphere, namely: 

X = Y^SP-^x&p) = e^Y^ST-^xi^p) . (7.17) 

Notice that we have assumed that x does not depend on the coordinates of the extra dimen- 
sions. Moreover, in (7.17) the functions Yi(S 7 ~ p ) are the scalar spherical harmonics on the 
S 7 ~ p sphere which, among other quantum numbers, depend on an integer I. The Yi(S 7 ~ p ) 
are eigenfunctions of the Laplacian operator on S 7 ~ p , with an eigenvalue which depends on 
/ and p: 

n S 7- P Y t = -l(l + 6-p)Yi. (7.18) 
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By plugging this ansatz for %, the action (7.14) becomes the action of a free massive particle 
propagating in AdS\ +q 

M 2 



q(x) 



-T D(7 + q _ p) 7V ; / dtd q 1 (dp ^/-g A ds 1+q 
where jVj is the integral 



Mi 



d 7 ~ p ( {Yi(S 7 ~ p )) 2 , 



and the mass M 2 of the I th Kaluza-Klein mode is given by 



M 



2 _ 1(1 + 6- p) _ 41(1 + 6 -p) 
u 2 (5 — p) 2 



(7.19) 



(7.20) 



(7.21) 



The equation of motion of x can be simply written as the equation of motion of a free massive 
particle in AdS q+ \ 

(-□ ff+1 + M 2 ) x = 0, (7.22) 

with being the d'Alembertian in AdSq+\. In Euclidean signature, the regular solution 

of (7.22) can be written in terms of the modified Bessel function of the second type 



X(t,p) = e iEt x(E,p) = e^(Epf/ 2 K & (Ep), a d = f \jM 2 + q - = 21 ^ p P , (7.23) 
and the asymptotic behaviour of x as the boundary is approached turns out to be 



X ~ P 



2a- (~{0) 



xL 0) + ...)+ P 2a + (xf + ...), 



(7.24) 



with the exponents q t obtained by solving 

- 2a (2a -9), 



M 2 



OL- 



I 

' 5-p' 



5~p 



(7.25) 



Notice that the parameter a in (7.23) is nothing but the difference between a + and a_ (i.e. 
a = a+ — a_). 

In the framework of holographic renormalisation, the non-normalisable fluctuation x^— 
acts as a source of the operator O x dual to the fluctuation Xi whereas x+ as a vev - 111 
this respect some comments are now in order. Sourcing the mode x m the bulk theory is 
equivalent to sourcing an irrelevant operator O x at the boundary. We are working pertur- 
batively in the mode Xi X being a fluctuation. This is equivalent to consider parametrically 
small sources of the irrelevant operator O x . As first noticed in [29] (see [49] for a further 
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reference on holographic renormalisation of irrelevant operators), one can consistently holo- 
graphic renormalise the theory order by order in x, making the n-point correlator of the dual 
irrelevant operator normalised up to some fixed n 5 . We are interested in computing the two 
point function, and therefore it will be enough to consider and renormalise the part of the 
action which is quadratic in the fluctuations. 



7.1.1 Holographic Renormalisation 

As a first step, let us consider the asymptotic expansion of x in a neighbourhood of the 
boundary p = 0. This expression can be written as 



- P -* 2a 

X ~ P 



k=0 

oo s 
k=0 1=2 



E~(2fc) 2k , r o(2fe) 2k 1 

X+ P +5/3,5-1 P lo g/° + 



.fe=0 



k=0 
P = 



5 — 1, a £ Z 
[a] otherwise 



(7.26) 



where [a] denotes the integer part of a and the logarithmic terms can be present if and 
only if a is an integer (implying that p must be even), which can never occur at leading 
order as it can be easily seen from the explicit expression for a in (7.23). More precisely, the 
logarithmic terms can start to appear at order 0(p 2a ~ +n ), with n being a (non-zero) multiple 
of 4. Therefore, the modes characterised by logarithmic terms in the asymptotic expansion 
close to the boundary have quantum number I given by 



I = (h-2p)n- (3-p), 



3 — p 
n > £z e 



p 



0, 1, 2, 



(7.27) 



5 - 2p 

where p and h are respectively related to p and n by n = 4n and p = 2p. Furthermore, 
in absence of the logarithmic terms, the coefficients in such an expansion can be obtained 
recursively to be 



~(2fc) 



□ g x| fc - 2) . r(d - fc) 



^□, fc x' 0) 



k G [1, +oo[, 



(7.28) 



Ak(a - k) 4 fc r(A: + l)r(fi) g AT ' 
with □ being the d'Alembertian with respect to the boundary metric g. In presence of the 
logarithmic terms, the recursive relation (7.28) still holds for the coefficients X- \ k < a, 
and the coefficient i?^ of the first logarithmic term is 

I — I ~(2S-2) r>i — ifi ~(0) 



2a 



4°r(d + l)r(a) 



The on-shell action for the fluctuation is 



5 M ~ \ [ dt<p-\ Jgxpd p 

^ J M F 



X 



(7.29) 



(7.30) 



3 We thank Bait van Rees for discussion about this point. 
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with e being a cut-off for small p and M. e is the boundary of AdS\ +q at p = e. Using the 
asymptotic expansion (7.26), the divergent part of the action turns out to be 



5<x) ^ l I dtd q ~ l (e 

Idiv 2 ' 



1/- _-2a 



i9 fc 

2 E e2fc Z («- + k - m ) *- m) *- fc " 2ra) + 

fc=0 m=0 



(7.31) 



+2(a + + a_)xL 0) ^ ) e 2 "log £ ], 



and by re-expressing the action in terms of the original field x, one gets the following coun- 
terterm action 



J Mr 



P 



k=0 



(7.32) 



with a (0) = 2a_. Thus, the renormalised action becomes 
5 (x) | = lim \S (X) \ +S^} = 

Iron g. ^.q L \£ cl j 



lim 



L D(7- <J + p) 



M [ dttf-^Jd. (xsd £ x)\ + 

J M-f V ' L 



+ 



2d(7 



(7-9 + p) 



-V, / dtd^\J^ e 

J Mr 



P 



Y J « i2k) X^X + a^xnfx\oge 

k=0 



(7.33) 



where A/"; has been defined as the integral of the spherical harmonics in (7.19) 



7.1.2 One-Point and Two-Point Correlation Function 

Let us now use the renormalised action to compute the one-point correlation function, which 
is defined as 



lim 



1 1 6S^\ 



2» 4 



The first variation of the renormalised action turns out to be 



6 S (x) \ 



D(7-q + p) 



and therefore the one-point function is given by 

(6 X ) = r^M [2(a--a+)xf + 8 p ,a-i0f] 

T (p, n) -A/j 



"2 : X+ + 0/3,5-1 



2d? X 



(0) 



5 — p 



4 5 r(a + l)r(d) 



(7.34) 



A/1 / dt^Ce- 2 "- [2(a--a+)xf+#f] 5x(t,e), (7.35) 



(7.36) 



Typically, the second term in (7.36) can be removed completely by adding a finite counter- 
term, corresponding to the matter conformal anomaly in AdS q +\ [50]. 
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The coefficient x+ can be read off from the full-solution (7.22). In order to perform 
this calculation, let us extract the time dependence of x(t, p) and work with x(E, p) (see eq. 
(7.17)). Expanding the latter in a neighbourhood of the boundary we have to distinguish 
two cases, namely 



a non-integer: 
x(E,p) = 

a integer: 

x(E,p) = 



„2»_ 



1 + 



H a ~ 2 7 e 



(Ep) 2 
4(5 + 1) 

(Ep) 2 
4(1 - a) 

2 



+ 



+ ...+ 



(EpY a Y{l-a) 
2 2 « T(l + a) 



+ 



2a 



(~l) a (Ep) 
2 2& T(a)T(a - 1) 



(7.37) 



i E 
l) l0g 2 



(-1) 5+1 2 {Ep) 2& , 

+ 2^r(d + i)r(d) log/5 + 



where 7 EM = -577 is the Euler-Mascheroni constant and, for a given integer n, H n is the 
harmonic number (H n = Ylk=i t)- By using these results the coefficient x+ turns out to be 



E 2& r(l-a) -,(0)/pn 
' 2^r(l+a)^- 



-(0) 



(-l) a .E 2 



2 2ci r(«)r(a-i) 



27-^logf ^(^) 



i(fi-l) 



a non-integer 
a integer 



(7.38) 



The correlator (7.36) therefore acquires the following form 



2l + 6-p 



4i+3(5-p) 



V 



where A/", (p ' 



n) 



} = -Ni P ' n) 2 Z — [g 2 a (E)]^E^-e(E) 

' gl M NE p - 3 and Q(E) 



(7.39) 



-M T( v ,n) (T( P ,n) has been defined in (7.9)), g cB (E) 
E~ 2a x ( + • Notice that Q(E) can be written in terms of the source X- (E) by using (7.38). 
This relation depends on whether a is integer or not. As a = 1 + ^fr^-, one has that a 
is integer iff 21 + 1 = (mod(5 — p)). This never occurs if p is odd and it happens in all 
cases if p = 4. In the remaining cases p = 0, 2 it only occurs for some particular values of the 
Kaluza-Klein mode I. Furthermore, for a integer, the only relevant term is the one containing 
log E, while the others are scheme dependent and we will omit them. 

Finally, the two-point correlator can be obtained from (7.39) by differentiating with re- 
spect to the source %_ . It is easy to see that for a non integer, it is just given by contact 
terms 

ioA) = -A?--'! i&m * **»*S^i> (7 . 40) 



r(i + a) 



while for a integer it has a logarithmic dependence 

(OA) 



2(-l) a A/"/ P ' n) 4 i+ 3(5- P ) £ 23 

Yy k ff (^)] ^ 5-p log 



4«r(a)r(a - 1) 



2a • 



(7.41) 
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The two-point correlator in position space can be obtained by Fourier transforming (7.40) 
and (7.41). With this purpose, we will use the relations 

I a (t) d = f J dEe-* tE E^ = 2 2fi + 1 ^ r( r ^ 2 _~ ) " ) |t|1 1 +25 , a non-integer, 

itBrfch™. jp _ ,_ ni+a ( 2q ! t, ( 1 



I b (t) = J dEe'^E^logE = , d integer, 

(7.42) 

where 7£ (•) indicates the renormalised version of its argument and is defined as 

n ( --^J] = rn ^ - p^-5^\t) , (7.43) 
/ |i| 1+2a (2d)! w v ' 

with 5^ (t) denoting the m th derivative of the Dirac 5-function. The two-point correlator in 
position space becomes 



<dA> = -A/^4dv ^ " ^ ^ " [y Zl , d non-integer, 



■r (1/2 - d) r (i - d) [g 2 ff ( t)] 
r(-d)r(i + d) |t| 

(7.44) 

= — ' , ' ; -72. , IOA , a integer, 

4°r(d)r(d - 1) I |t| 2A * 

where the effective coupling constant is defined as g 2 ff (i) = g 2 M iV|t| 3-?) and A/"/ p ' n) ~ N. 
Notice that the exponent of \t\ in (7.44) provides twice the generalised conformal dimension 
A x of the operator O x , which is given by 

A v = — I + 2 . (7.45) 

For p = 3 this result coincides with the correct value A x =1 + 2 for the truly conformal 
case [20,24,25]. Notice that A x is fractional for p < 3. 

7.2 Fluctuations of the Angular Embedding Function and World-volume Gauge 
Field 

Let us move on to the angular embedding, whose fluctuations are coupled to the ones of the 
world- volume gauge field. We analyse the fluctuations around the configuration 6 = #( p , n )j 
namely 

8 = 0<p,») +£, = ^ + -^S , (7-46) 

where the only non-zero component of J-^ is = cos0( Pin) and, in terms of the rescaled 
dilaton (p defined in (7.6), the fluctuation field can be rewritten as 

= fafi . (7.47) 
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When these fluctuations are switched on the induced metric on the probe D(8 — p)-branes 
acquires a further term given by 



~(i) 2 

9Wti = u v 



{d a £){dp£) + (2 sin cos6> (:p , n) £ + (cos 2 6» {p>n) - sin 2 6» {p>n) ) ) £ 2 6 a a 6p g ab 



(7.48) 

where a, b are indices along the S 7 ~ p sphere and g ab denotes its round metric. The D(8 — p)- 
brane action for these fluctuations acquires the following form 



p— 3 

5 D(8 _ P) = -T (p>n) I rf 9_p (e^y-(/ 4li;;2 ^ S T-,x 

r (0 

Vdet{I + X}+ V ~ tp 



{•Ftp + Ftp) Cp(^(p,n) + 



sin 8 p 6 {p , n) 



(7.49) 



where C p (9) is the function defined in (3.18) and X is a matrix defined as 



(7.50) 



Expanding the action (7.49) and keeping the terms up to the second order in the fluctuations, 
one gets 



T(p, n) 
,2 



d 9 P C e'-pVyZ-g^y/g^X 



(P, n) 



3-p, 



/a/3/75 



7-p Zes-p* f t p 
V~9Ads 2 sin 3 

17 (p, n) 

(7.51) ' 



where the open-string metric t7 Qi a is defined as in (7.3). Such a Lagrangian coincides with 
eq (2.61) in [20] if the components f ab along the S 7 ~ p directions are taken to be zero. Let 
us represent f a p in terms of a potential a a as f a p = d a ap — dpa a . Then, the equations of 
motion for the fluctuations (£, ap) have the following form 

3-P, 











1 



=d a V-9G ap dp^ 



-9 



+ 



7-p 



Up sin 2 #( Pj „j 



e 5 '" / 



tp 



\f~~9~Adsl Sm "(p, n) 



1 



8111 



(P, ") 



(7.52) 



. p — 3 

where \/— 5 = e 5 -? ^ ' yf '— g AdS2 ^/g s 7 -p ■ Similarly to the case analysed in Section 7.1, the 
equations of motion (7.52) are equivalent to the ones in AdSq + \ compactified on a T q ~ l and 
with q analytically continued to take the value 

2 



5 — p 



(7.53) 
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Indeed, it is straightforward to check that one can rewrite the system (7.52) as 

ftp 



yf 9Ads i+1 srn "(p, n) 







1 



y/ 9AdS i+1 
u- 2 



0: 



7 — p 



sm( 



*X - sfsf 



+ (7.54) 



V9^9 ac G BD e 2 ^f cD , 
Let us now analyse in details the different modes, following the classification in [51]. 



7.2.1 Coupled Modes 

Let us now choose the following ansatz for the scalar £ and the gauge field 

i = e iE *y z (5 7 ^)|(p), a a = e^VaY^-PHp), a p = e^Y^^a^p), (7.55) 

with all the other components of the gauge field set to zero. Notice that the field component 
a a is a gradient in the S 7 ~ p sphere. Therefore, it is always possible to perform a gauge 
transformation such that the components of the gauge field along the S 7 ~ p are set to zero 
while the non-zero components are the ones along the AdSVdirections. Accordingly, in the 
following we will adopt the following ansatz for the components of the gauge field 



at = -iEe^YiiS^atip), a f 



JEt 



(7.56) 



With such an ansatz, the equations of motion (7.54) acquire the following form 

n .„ 







(5-p) 
.7 — p 



= E 2 



= E 2 



y/-9Ads i+1 sml 



(5-p) 2 
41(1 + 6 -p) 



y/~9AdS~^ Srn 9(p, n) 

9p P at 



E 

.7 — p - 



ts 4+1 sinl 



^9a 

9 PP at 



\J~9Ad 



E 
+ d p 



+ 



(5-p) 2 ^/-9 A ds i+1 sin {Pi n) 
41(1 + 6 -p) gtt(a p -d p a t ) 



(7.57) 



(5-p) 2 y/-9Ads s+1 sm0 (p , n) 
9tt (a p - d p a t ) 



y/~9AdS q+1 s i n 6(p, n) 



where + ! is the d'Alembertian operator for AdSq+i. This is a system of four equations in 
three unknowns at, , and therefore one of the equations must be redundant. It is easy 
to see that solving the second and third equations in (7.57) for their last terms and inserting 
such solutions in the last equation, the latter turns out to be an identity. So, we can consider 
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the last equation in (7.57) as redundant and focus on the first three. Their form suggests 
that it is convenient to define two new fields ( £, f] 



1= ~ie iEt 



.7 — p i 



£ \e iEt . 



(7.58) 



\ y/-9Ads 4+1 sine £ 

Now, eliminating the field at by using the second equation in (7.57), the equations of motion 
can be reduced to the following system of two equations 



(5-p) 



1_ [/(/, + 6 - p) + (7 - p)(6 - p)] | + t?, 



(5-p) 5 



4(Z + 6-p) 4(7 -p)Z(Z + 6 -p) 

U ~~ U 'q + l'l te _\9 / ' lr __\9 V 



(7.59) 



(5 — p) 2 (5 — p) 2 

It is convenient to rewrite the above system of two equations in matrix form 



n 9+1 tp = M(i, p) tp, v 



(7.60) 



and the mass matrix _M( iiP) , as in [20], is 

fl(l + 6-p) + (7 -p)(6-p) p-7 \ 



(5 - pf 



(7.61) 



y (p - 7)1(1 + 6 - p) 1(1 + 6 -p) J 

The eigenvalues A (l) and eigenvectors ip^ of M.(i,p) are 
, 1} _ 4(/ + 6-p)(/ + 7-p) 



A 1 

A (2) 



(5-p)2 
4/(/ - 1) 



(5-p) 2 ' 

and they satisfy the differential equation 

D i + 1 iJ) w = A (i V (i) - 



= (/ + 6 _p)£_^ / > o, 
^ (2) = | + ?7, I > 1, 



(7.62) 



(7.63) 



It is now straightforward to obtain the behaviour of the normalizable and non-normalizable 
modes for tp^ in a neighbourhood of the boundary p — > 0. One gets 



( (!) 



Z+6-p 
5-p ' 



Q: 



(1) . 


l+7-p 


+ - 


5— p 


(2) = 


l-l 




5-p 


(2) _ 


I 


+ — 


5— p ' 



(7.64) 
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In order to compute the correlation functions for the decoupled modes one can notice 
that their equations of motion (7.63) are just the equation of motion for free massive scalars 
in AdSq+i and, therefore, the dynamics of the fluctuations - up to quadratic order - can be 
described by the following effective Euclidean action 



(7 + 8 -p) 



1 



(7.65) 



The holographic renormalisation of this action goes exactly as in 7.1.1, so that the one-point 
function can be written as 



Pi 

def (i) (i) 



(7.66) 



with d w = a\_' — a_ and Q(E) given by 



G(E) 



r(l-qW) 

r(i+aW) v - ' 



a [l> non-integer 



(-1)' 



k 4 aW r(aW)r(aW-i) 



H 



27e 



2 i E 



d w integer. 



(7.67) 

From the values of a± written in (7.64), one immediately determines the values of d w , 
namely 



21 + 3 



+ 2 , 



d (2) 



21 - 1 



(7.68) 



5 — p 5 — p 

Then, it follows that d (1) is integer iff 21 + 3 = (mod(5 — p)), while d (1) £ Z only when 
2/ — 1 = (mod(5 — p)). As in the case of the \ fluctuation, this can only happen when p is 
even and always happens if p = 4, whereas for p = 0, 2 it occurs for some values of L 

Similarly to the case of the fluctuations of the Cartesian coordinates, the two point func- 
tion turns out to have the following functional form 



[g 2 eS (E)] «-* E 



[^^^^lOg., 



d w non-integer, 
q w integer, 



(7.69) 



which in the coordinate space becomes 

I A \ [#eff(*)] 



1 

5-p 




d w non-integer 



d w integer. 



(7.70) 
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with the exponent A ( ,*' given by 



A ( , l) 



Q W + 1 



5 — p 



2 (J + l) + 3, 



5— p 



K 5-p 



2 a - 1) + 1, 



2. 



(7.71) 



Notice that the exponents correctly reproduce the conformal dimensions in the case 
p = 3, where A^ =1 + 4 and A^' = I [20,24,25]. Again the conformal dimensions are 
fractional for p < 3. 



7.2.2 Internal Gauge Field Modes 

Finally, let us consider an ansatz for which the scalar £ is set to zero and the only non-zero 
component for the gauge field are the ones along the S' 7_p -directions 



0. 



fta = d t a a , fpa = d p a a , fab = d a ab — dba a , 



(7.72) 



with 



a a = Y a {S 7 ~ p ) a(t,p) ee Y* {S 7 ~ p ) e* Et a(p), (7.73) 

where Y a (»S 7_P ) is a vector spherical harmonic on the S 7 ~ p sphere. With such an ansatz, 
the equations of motion (7.54) can be written as 

= d t {V a a a ) , 
= d p {V a a a ) , 



9 PP 



V 9Ads i+1 
9tt9 P p 



dfa b - a 



9tt 



V 9Ads i+1 
u~ 2 (Aa b - V 6 V r a c 



aV 



(7.74) 



"p 

where the operators A (the Hodge-de Rham operator for one-forms on the S 7 ~ p ) and V& are 



defined as 



Aa 6 



(7.75) 



with being the Ricci tensor on the S' 7_p -sphere. From the first two equations (7.74), one 
can deduce the requirement 



V h a b 



0. 



(7.76) 



which implies that the last term in the third equation (7.74) vanishes. Furthermore, from 
the decomposition (7.73) the condition (7.76) implies that the vector spherical harmonic 
Y a (S 7 - p ) has to be such that 

V b Yl (S 7 ~ p ) = 0. (7.77) 
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Considering also that the vector harmonic (5 7 p ) is an eigenfunction of the Hodge-de 
Rham operator A defined in (7.75) 

A Yj (5 7 ^) = -(I + + 5 - p)Y l a {S^p) , (7.78) 

the third equation in (7.74) takes the form 

= (a+i - Ml) a, (7.79) 



where D s +i is the Dalambertian operator in AdS s +i with s = 2 — q = 2(4 — p)/(5 — p) and 
the mass M~ is 

M l = 7T^2 [(l + m + 5-p)]- (7.80) 



a 



The behaviour of the normalizable and non-normalizable modes a ~ a_ p 2a - + a+ /) 2a + in a 
neighbourhood of the boundary is given by 



Mf = 2a(2a-s) 



«- = -5^( Z + 1 ) ) 

(7.81) 

a + = 5h( l + 5 ~P) 



In this case 

21 + 1 

a a = a + — a_ = + 1 . (7.82) 

5 — p 

Proceeding in a similar fashion to the analysis in the previous cases, it is easy to find that 
the two-point function of the operator dual to these fluctuations has the following form 

i 

5— p 



I A A \ [fleff(^)] 5 

{O a Oa) ~ — 1^|2^ — , a a non- integer 



(7.83) 



with A a being given by 



A a = 5 a + l-_ L- = -J— l + 1 (7.84) 
5 — p 5 — p 



Notice that A a = A x and for p = 3 the exponent A a coincides with the conformal dimension 
of the operator O a in the AdS§ x 5 5 background, which is A a = 1 + 2, as it should [20,24,25]. 



8 Conclusion 

In this paper we investigated the insertion of impurities in (p+l)-dimensional Supersymmetric 
Yang-Mills theories, which are (trivially) non-conformal. The ambient theory is holograph- 
ically described by the near- horizon geometry generated by a stack of iV Dp-branes, while 
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the impurities are added by introducing probe D(8 — p)-branes in this background in such 
a way that the induced world- volume metric is conformally AdS2 x S 7 ~ p . The background 
RR (7 — p)-form potential induces an electric gauge-field on the world- volume of the probe 
branes, giving rise to a bundle of strings stretching in the radial direction and forming a flux 
tube. 

We analysed in some detail two possible classes of configurations for such systems. In 
the first one, the flux tube is straight - which corresponds to a constant embedding function 
for the probe branes -, while in the second one two flux tubes get connected in the bulk. 
From the gauge theory point of view such configurations respectively correspond to single 
impurities and to dimers. 

For these systems we studied the basic thermodynamic properties such as the free energy, 
entropy and specific heat. Interestingly, the impurity entropy turns out to be in general 
non-analytic in the filling fraction, except for the case of the D4/D4 systems for which we 
had been able to find a closed form. This latter system is actually the only one showing a 
positive specific heat: for p = 3, the specific heat is zero, while for p < 3 it is negative. This 
is a signature of a thermodynamic instability. We also computed the impurity susceptibility, 
which is constant for p = 4. 

We then analysed the case of the hanging flux tubes, for which the position of the probe 
branes is fixed in the transverse direction, while their embedding in the conformal-AdS' man- 
ifold is controlled by a scalar. At finite temperature, there are two possible configurations: 
in one the two flux tubes are connected in the bulk and lie outside the black hole, while in 
the second one the two flux tubes end into the black hole. The transition from the second 
configuration to the first one is of first order and it corresponds to the dimerisation transition. 
We analysed the thermodynamics of these connected configurations. All the thermodynamic 
functions (free energy, entropy, internal energy, specific heat and latent heat) can be written 
in terms of hypergeometric functions. We studied the competition between the two con- 
figurations and determined the temperatures at which the dimerisation occurs. This phase 
transition turns out to be of first order. Furthermore, the specific heat for the dimer con- 
figuration is always positive for any p < 3 and vanishes as T — > 0. As mentioned before, 
the specific heat for the straight flux-tube configurations is negative for p < 3. One possible 
interpretation of this result can be that such a configuration is just not allowed for p < 3 
and that the flux-tubes are forced to reconnect and form a dimer. 

We further studied the stability of the systems at zero temperature by analysing the 
fluctuations of the probe branes in the straight configuration (impurity fluctuations). The 
fluctuations decouple into two channels: one contains just the fluctuations of the Cartesian 
coordinates, the other channel instead contains the coupled fluctuations of the angular em- 
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bedding function and of the world-volume gauge field. The modes in the second channel can 
be conveniently decoupled. Interestingly, all these modes satisfy the equations of motion of 
free massive scalars in a higher-dimensional AdS-sp&ce, where the enhancement of dimen- 
sions is due to the presence of a dilaton with a non-trivial profile. These modes are dual to 
irrelevant operators. In order to make sense of the correlators of such operators, one needs 
to holographic renormalise the action for these modes perturbatively up to order n, if one 
needs to compute the n-point correlators. For the aim of characterising the theory, it may be 
enough just to compute the one- and two-point correlators. Therefore we needed to renor- 
malise the quadratic action which, as we have mentioned, is just the action for a free massive 
scalar in a higher-dimensional AdS-sp&ce. Then, we computed the correlators and wrote 
them in such a way that the underlying generalised conformal symmetry is manifest. From 
such an expression we can extract the generalised scaling dimension for all the operators we 
are interested in. Remarkably, such generalised scaling dimensions are fractional for p < 3 
since the term which depends on the KK quantum number I gets multiplied by the factor 
2/(5 — p), which is just 1 for p = 3. It would be interesting to understand this scaling from 
a purely field theoretical point of view. From the discussion about the generalised conformal 
structure, the most naive expectation is that the impurity action has a factor with a suitable 
power of the dimensionful coupling constant. 

It would be interesting to repeat this fluctuation analysis for the case of the dimer con- 
figuration. One would expect that some of the dimer fluctuations satisfy the Heun equation, 
as it happens in the ABJM case [26], which is connected to integrable models. If on one side 
the emergence of this integrable structure could be understandable for the case p = 3, it 
might be a bit counter-intuitive in the general p < 5 system. However, as we discussed, the 
modes for p ^ 3 can be generally seen as modes propagating in a higher-dimensional AdS, 
so that the features of these systems may be deduced from conformal systems. 

As we mentioned in this paper, while it has been showed [47] that the introduction of D5- 
branes in a D3-brane background such that the world- volume induced metric is AdS2 x S 4 is 
dual to a 1/2-BPS Wilson loop operator in N = 4 SYM in the anti-symmetric representation 
of the gauge group U (N), it is not yet clear whether a similar interpretation holds in the case 
of D(8 — p)-branes in the Dp-brane background with the same type of embedding. Therefore, 
a natural thing to do would be to clarify this issue. In this spirit, one can also study fermionic 
excitations in the Dp/D(8 — p)-system as done in [25] for the conformal case p = 3. 

The analysis we carried out in this paper can be generalised to the case in which the probe 
D(8 — p)-branes wrap different submanifold of S^ 8 ~ p \ as in [52]. More interestingly, one can 
think to introduce chemical potential to the system by considering Reissner-Nordstrdm type 
of backgrounds and analyse the configuration allowed and the phase structure. Importantly, 
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one could compute the resistivity as a function of the temperature and to study whether the 
holographic impurities induce a minimum similar to the one that occurs in the condensed 
matter models reviewed in section 2. 

A final - and longer term - direction concerns the effect of the backreaction of the D(8— p)- 
branes. Presumably, this would require dealing with a bubbling geometry, generalizing the 
one found in [24] for the conformal p = 3 case. 
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A Polar angles for stable embeddings 

In this appendix we analyze the main properties of the polar angles of the flux tube configu- 
rations. We begin by listing the different functions A P)7l (0) for < p < 5, namely 

2r « / ~ . a n . . o « . ~ \ ~ / ~ 7i 
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As mentioned in the main text the functions A Pin (6) depend on the quantization number n 
and on the number of colors through their ratio v = n/N (the filling fraction). Moreover, 
one can check explicitly from (A.l) that these functions satisfy 



A p ,n(0) = -A Pi 7v_ n (7r - 6) , 



(A.2) 
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from which (4.5) follows immediately. As shown in [20], the A p ^ n (9) are monotonically in- 
creasing functions of 9 in the interval < 9 < tt. Moreover, one can check that: 

WO) = 4^ £ , A P A*) = 2^ ±±± ( 1 - £ ) • (A.3) 

and, therefore, A Pjn (0) < if n > and A p .„(7r) > if n < N. Thus, it follows that there 
exists only one solution 9( p , n ) £ (Oi 71 ") °f the equation K Pjn {6(r, jV \) = for each n in the 
interval < n < N, i.e. there are exactly N — 1 angles which correspond to non-singular 
wrappings of the D(8-p)-brane probe on the S 7 ~ p sphere. Let us work out the expressions of 
some of them for different values of p. In the conformal case p = 3, the angles 9^ n ^ are the 
solutions of the equation: 

— — — Tl 

0(3,n) - cos 0(3,n) sin d (3,n) = K ■ (A.4) 

For p = 4 the angles 9u n \ can be immediately obtained from (A.l), with the result 

— IX 

cos# (4>n) = 1 - 2— . (A.5) 

The analytic expression of the angles 9( Ptn ) for p < 2 is more difficult to obtain. Let us 
analyse in detail the equation which determines the angles for p = 2, namely the solutions 
of the equation A2 n (6) = 0. In terms of the filling fraction u, after using the third expression 
in (A.l), this equation becomes 

cos 3 £( 2i „) - 3cos0( 2 , n ) = 4^ - 2 , (A.6) 

which is a cubic equation in cos#( 2 ,n)- I* 1 general, an equation of the type 

x 3 + rx = s , (A. 7) 

can be solved for x by means of the so-called Vieta's substitution, namely: 

x = w — . (A.8) 

3w 

Indeed, by substituting (A.8)) into the cubic equation (A. 7), one obtains the following 
quadratic equation for u; 3 : 

(^ 3 ) 2 " Hv 3 ) - ^ = 0. (A.9) 
In our case x = cos 9i 2 n ), with r = —3 and s = Av — 2. The two solutions for w 3 are just 



3 - 2v- 1 ±2iy/v(\-v) . (A. 10) 
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Remarkably, the right-hand side of (A. 10) is a complex number of modulus one and, therefore, 
can be represented as 



w 3 = e ia 



w = e a , 



where a is an angle such that 

cos a = 2u — 1 , 



sin a = ±2 y z/(l — i/) . 



(A.ll) 



(A.12) 



Notice that the relation between x = cos 6 and w is, in our case, given by (A. 8) with r = —3, 
namely 



It follows that 0(2,n) is given by 



COS 9( 2 ,n) = 2 cos ( — 



X = w + w 
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37 ' JV 
As a further example, let us point out that the angles for p = 1 can be obtained by solving 
the following transcendental equation 

n 



N 



7T . 



(A.15) 



B Dimer integrals 



The purpose of this appendix is to derive the integrals needed in the calculations of section 
5 of the thermodynamic properties of the dimer configurations. First of all, we define the 
following two integrals Ii(a,7) and /2(«i7) as 

roo j 

T , \ del I UZ 

Il(a,7) = 



h (a, 7) 



V0 a -l)(* a -7) ' 
z a dz 



(B.l) 



'1 yj{z a - l)(z a -7) 3 

where a and 7 are real numbers such that a > 2 and | 7 | < 1. These integrals can be 

performed in terms of hypergeometric functions, namely: 

1/ 1 1\ /l 13 1 
Ua,^) = -B[l ,- )F -, 1 ; 57 



1 / 1 1\ /3 13 1 
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(B.2) 



where B[x,y) is the Euler Gamma function : B[x,yj = T(x)T(y)/T(x + y). Let us next 
consider the integral: 

-R 



J(a, 7) = lim 



zr — 7 
z a - 1 



dz - R 



(B.3) 
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To evaluate this integral we shall proceed as follows 6 . First of all, let us rewrite the integrand 
in (B.3) as 



7 d 



z a - 1 dz 



z a - 1 
z a — 7 



a(l -7) 



2 ^(z" - l)(z° - 7 ) 3 
1-7 



(B.4) 



y/(z*-l)(z°< -7) 

Using this result it is straightforward to relate J(a, 7) to the following combination of the 
integrals Ii(a,7) and /2(a) 7) defined above 



a 



J(a,7) = (1-7) h(a,j) - h{a,i) 



(B.5) 



Moreover, by using identities satisfied by the hyper geometric functions, one can show that 
(for a > 2) 

, (a , 7) = -^ B ( 1 -Il) F (l, _I 1 1). (B.6) 
2a V a 2/ \ 2 a 2 a / 
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